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ABSTRACT 


A possible measurement program designed to obtain 
the information requisite to determining the feasibility 
of airborne and/or satellite-borne LDV (Laser Doppler 
Velocimeter) systems is discussed. Measurements made 
from the ground are favored over an airborne measurement 
as far as for the purpose of determining feasibility is 
concerned . 

The expected signal strengths for scattering at 
various altitude and elevation angles are examined; it 
appears that both molecular absorption and ambient tur- 
bulence degrade the signal at low elevation angles and 
effectively constrain a ground based measurement to ele- 
vation angles exceeding a critical value. 

The nature of the wind shear and turbulence to be 
expected are treated from a linear hydrodynamic model - 
a mountain lee wave model. The spatial and temporal cor- 
relation distances establish requirements on the range 
resolution, the maximum detectable range and the allow- 
able integration time. 

The ability of available LDV instrumentation to 
measure typical and shear profiles and turbulence levels 
is then assessed. 

Analysis of the operation of LDV systems under 
very general uontli uions is given in appendix I and a 
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mechanical shutter described in appendix II is proposed 
for solving the ringing problem for the pulsed LDV system. 
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I . SUMMARY 

There are presently under development coherent laser 

detection systems designed to measure atmospheric winds 

(121 

and turbulence at both short and long ranges , A 

continuous wave (CW) CO 2 laser system has been demonstrated 

to provide v?ind shear infomation at altitudes up to 500 feet 

(2 3) 

and ranges to 1000 feet . A pulsed version of tlae same 

instrument designed to detect clear air turbulence (CAT) has 

been mounted on an airborne platform and operated at various 

(4) 

altitudes up to 40,000 feet . Clear air detections at 
ranges up to several miles have been obtained at lower 
altitudes (less than 15,000 feet) . Clear air returns were 
not consistently obtained at higher altitudes and it v;as 
concluded that greater sensitivity or more output power 
would be required for operation at these altitudes. 

The objective of the latter program to develop a 
system capable of detecting hazardous regions of turbulence 
(CAT) for installation in commercial aircraft, v/hareas the 
former program was oriented toward the development of 
scanning velocimeters that v/ould provide data on local v/ind 
shear for various application including aircraft operations in 
airport environment, artillery applications, etc. In addition 
to these objectives, there are a number of other important 
applications of such remote wind shear measuring systems. 


In particular, the ability to measure upper atmospheric 


winds on global scales from a satellite platform is sufficiently 
attractive for weather forecasting and meteorological studies 
to merit quantitative analysis of the potential of such a 
measurement. At a more local level these systems can be 
used to measure the v/ind shear field in and near various 
meteorological phenomena (fronts, large amplitude internal 
waves, large scale convection cells, regions of developing 
storms, etc.) . 

Preliminary tests with a pulsed coherent airborne LDV 
(laser doppler velociraeter) have been inconclusive, and it 
has not been possible to assess the feasibility of measuring 
winds and wind shears of long ranges from these tests alone. 
Uncertainties exist both in the level of the aerosol back- 
scatter at high altitudes, as well as in the instrument per- 
formance. Attempts to sample the aerosol content using 
mechanical samplers have been carried out also v-?ith incon- 
clusive results. 

In this report we will discuss a possible measurement 
program designed to obtain the information requisite to 
determining the feasibility of airborne and/or satellite- 
boriie laser doppler systems. We will argue that for the 
purpose of determining feasibility and for the acquisition 
of the requisite design information, measurements made from 
the ground are to be preferred over an airborne measurement. 




The advantages of a ground based measurement as 
compared to an airborne measurement are several and include: 

1) Availability of larger useful apertures (diameter of 
2 feet or greater as compared to a probable 1-ft. 
aperture for an airborne system) resulting in an 
enhancement of S/N at a given range of about 6 db. 

2) Availability of a relatively controlled laboratory 
environment for optics-laser-detectors and electronics 

3) Space and power limitations are not severe. 

4) The line of sight is stable. 

5) Ability to take extended and repeated measurements 
under a variety of conditions. Experiments can be 
repeated at relatively low cost. 

6) Ability to make extensive use of pulse integration 
techniques. 

7) Ability to utilize bread board laser and electronic 
instrumentation . 

The disadvantages of such a measurement include: 

1) Increased molecular absorption along the line of sight 

2) Degraded propagation characteristics by turbulence 
in the lower atmosphere. 

3) Poorer aspect for high elevation angles for measure- 
ments of the horizontal wind vector. 
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4) Inability to resolve fine scale vertical stratifi- 
cations. 

5) A lessened ability to identify turbulence from the 
spectral broadening of the scattered signal (greater 
interference from vertical gradients in the hori- 
zontal wind) . 

In the following sections we first examine (Section II) 
the expected signal strengths for scattering at various 
altitude and elevation angles. It will appear that both 
molecular absorption and ambient turbulence degrade the 
signal at low elevation angles and effectively constrain the 
possible measurements to elevation angles exceeding a critical 
value. This value depends upon output power, integration or 
averaging time and, at least for turbulence effects, ambient 
conditions , 

Follov/ing this analysis we treat in some detail 
(Section III) the nature of the wind shear and turbulence to 
be expected. The nature of the vertical profile and the 
spatial and temporal correlation distances establish require- 
ments on the range resolution, the maximum detectable range 
and the allowable integration time. 

The ability of available LDV instrumentation to measure 
typical v/ind shear profiles and turbulence levels is then 
assessed (Section IV) . Both ground and airborne platforms 
are considered. For a ground based measurement, it is 
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concluded that extensive pulse integration techniques 
may be needed to derive wind profiles up to 10 km, de- 
pending on the actual value of the backscatter coefficient. 
For detection of winds at the tropopause, the elevation 
angle is probably linxited to being greater than 30° . This 
implies a minimum vertical resolution of the order of 300 
meters . 

The current available pulsed LDV system^ ^has a 
ringing problem that limits the minimum detectable relative 
velocity. Since expected relative velocities at lower 
altitudes are relatively small, techniques for eliminating 
this ringing are required. One such technique using a 
mechanical shutter is described in Appendix II. 

In Appendix I, an analysis of the operation of LDV 
systems under very general conditions is given. Most of 
this material has been taken from a previous limited-dis- 
tribution report (Thomson and Dorian, 1967) , which is pre- 
sently unavailable and is included here for completeness - 
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II. ANALYSIS OF A GROUND-BASED LDV SYSTEM 


In this section the expected signal to noise ratio 
(SNR) for a ground-based pulsed LDV is evaluated. Expres- 
sions for the SNR derived from the general analysis of LDV 
operation in Appendix I are given. Atmospheric losses due 
to molecular absorption and turbulence are significant 
for the proposed ground based system even under good see- 
ing conditions, and limit the measurement ot elevation 
angles greater than 30 degrees. The sin.gle-pulse SNR's are 
small for the current system (pulse energy 20 m j , aperture 
diameter 30 cm); increased power (200 m j ) or larger optics 
(— 50 cm) V70uld be required to achieve a marginally detec- 
table signal on a single-pulse basis. Further improvement 
requires pulse integration. The ground-based measurement 
is found to have an SNR comparable to or slightly higher than 
an equivalent airborne system for similar conditions. A 
comparison of the ground-based LDV system with the Wallops 
Island radars, which have detected CAT produced by high- 
altitude internal v/aves, shows that the SNR's can be of the 
same order. 

A. S/N Expression s 

For a pulsed system operating in a collimated mode 
(focused at infinity) , the SNR is 


S/N = n 


J 

32 'hc/A' 


2 - 

D ng CT 



1 


where 


J - net output energy in the observation time 
D = aperture diameter 
T = pulse time 

no = particle density times scattering cross-section 
(mean value at range R) 

n = detection efficiency 

R = range 

L = propagation loss factor 


The factor L is the product of the molecular loss and 

the turbulent loss 

2 


At ranges less than the Rayleigh distance (R^ = irD /4\) 
the signal to noise ratio varies with range only as a result 
of the loss factor L; 


— 2 

S/N— n - ^ — L(R) for R < R* 

2if‘^hv D'^ 

For a 30 cm effective aperture diameter and \ = 10.6]J/ this 
expression is appropriate for ranges less than about 7 kra. 

Within the range limits, some benefit can be obtained 
by focusing on the range interval of primary interest- For 
pulse lengths less than the focal depth < 4R^A/D^^ , the 

signal to noise from the region near the focal point is 




s/N = n 


JX D no cT 


32hc 


L(R) 


R 


and exhibits an inverse square dependence on range in addition 
to the range dependence through the propagation loss factor. 

B. Loss Mechanism 

Molecular absorption results principally from CO 2 lines 
and the wings of the nearby water vapor bands. Two-way 
absorption losses as a function of the altitude at which the 
light is scattered are given in Figure 1 for various elevation 
angles of the line of sight. For a sea level installation, 
the losses are moderate (less than 12 db)' only for elevation 
angles (above the horizon) greater than 45°, An 
improved situation exists for a higher altitude installation 
under low humidity conditions. In Figure 2 two-way losses 
are shown for a source located at an elevation comparable to 
that at Boulder, Colorado (->' 5,000 feet). 

Under conditions of strong surface heating, low altitude 
turbulence will materially degrade the returns from the higher 
altitudes. The signal to noise ratio in this case is con- 
veniently expressed in terms of an effective coherent aper- 
ture radius r_ (r^ is related to Fried's effective aperture 

a a, 

8 

diameter r according to r = 3.18 r ) . When r„ is small 
o o a o 

compared to the actual aperture D, the optimum signal to noise 
(that obtained when the beam is focused on the range point of 
interest) is given by 


O 








In Figures 3 and 4 we have reproduced Fried's figures 
for effective apertures for both horizontal and vertical 
paths. These data may be scaled to finite elevation angles 
by the relation 

r^ = 0.0581 r3/5) 

where R is the slant range. For a vertically stratified 
atmosphere/ the following expression should be used: 



o 


where 0 is the angle from the vertical and H the scattering 
height. 

C. System Application 

The preceding relations have been used to derive 

estimates for expected signal to noise ratios for a suggested 

ground measurement. These values are shown in Table I- Under 

2 i 2 1/3 

moderate to good seeing conditions ( J dz < 7:<10 meters 

o 

turbulence is not a significant degrading factor for 30' cm optics 
In Table I we have assumed an overall optics-detector-electronics 
efficiency n of 2.0 percent. Theoretically, it should be 





TABLE I; Ground-Based LDV System 


Aperture diameter D = 50 cm 


^ -2/3 12 , „ nrL-12 1/3 

< 3x10 m ^ ; ICj^ dz < 7x10 m ' 


Detection Altitude: 10 km- 


Angle from Vertical: 45^ 


Pulse; 2\i sec (8u sec) 


na = 10 ® km ^ ster ^ 


Absorption = 6 db 


System/Detector Efficiency = 0-02 


Single Pulse SNR 

Pulse Energy 

S/N 

20 mj 

0.24 (0.8)* 

200 mj 

2.4 (8) 

Integrate 600 pulses 

incoherently 

Pulse Energy 

S/N 

20 mj 

6.4 (19)* 

200 mj 

64 (190)* 


Values in parentheses correspond to the sec pulse lengtl; 




possible in a pulsed system to achieve efficiencies approaching 
the detector quantum efficiency ( 30-50%) . Hov/ever, current 
estimates for the airborne system are much loiter and we 
choose 2% as a currently achievable value. Tv/o pulse lengths 
are considered; 2y sec and 8y sec. At an elevation angle 
of 45°, these imply a vertical resolution of 300 and 900 
meters respectively. For these parameters the single pulse 
SNR's are below detectability for a 20 mj pulse and are 
marginal for a 200 mj pulse. Since the external losses in 
this case are only 6 db, similar numbers (actually, a factor 
4 higher) will apply to the airborne platform measurement 
at the same range (14 km) . 

Under these conditions the single pulse signal to noise 
ratios are very low for the current system capability 
(J < 20 mj) and substantial increases in pulse output at 
power or sensitivity would be required for detection. Inte- 
gration (incoherent) of 600 or more pulses is required to bring 
the SNR up to acceptable values. It should be noted that 
these values apply to the entire signal and implicitly 
assume that the doppler spread in the returning signal is 
no greater than the inverse of the pulse width, that is, 

100 to 500 K Hz (or 0.5 to 2.5 meters per second). Substantial 
velocity gradients along the line of sight or significant 
turbulence levels v/ill broaden the signal beyond these limits, 
further reducing the SNR. 
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The above estimates imply that substantial pulse 
integration may be required to obtain useful data. The 
maximum allowable integration time is limited by the air 
transit time through the resolution volume. For .a ground- 
based system this time is 


I 



I 


.. cx sin6 

= T 

W 

where V is the horizontal wind velocity. For our nominal 
system with x = 5y sec and a wind velocity of 40 meters/sec, 
integration times of the order of 13 seconds are permitted. 
In the equivalent airborne system (travelling at a relative 
velocity of 300 m/sec) the integration time would be about 

2.5 seconds. At a pulse repetition rate of 200 pps and 

20 mj pulses (mean power of 5 KW ) , the ground-based system 
could average over 2600 pulses, yielding an SNR of about 

7.6 9 db) as compared to the airborne system value of 6.6 
Thus, we conclude that, vrhere operated under similar 

conditions, the ground-based measurement would achieve com- 
parable or slightly higher S/N than the airborne measurement 
The available increase in integration time in the ground 
system compensates for the increased absorption in the lov/er 
atmosphere. However, the signal to noise ratios predicted 
in this somev/hat conservative estimate are not high (less 
than 10 db) . For these estimates v/e have assumed that a 


I 
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spatial resolving power comparable to , 3 km to 1 km. 
was required. Improved SNR could be attained in the 
ground system by increasing the integration time to times 
of the order of the natural internal wave period of the 
stratosphere (typically several minutes) . Also, the pulse 
duration could be chosen to match more closely the velocity 
spread expected. The major source of CAT at high altitudes 
(see Section III) is presently surmised to be the instability 
of large amplitude internal waves near the tropopause. 

These waves and their instabilities have been detected 
by high power sensitive radar. In Table II we present 
comparative data for the Wallops Island radars and for the 
laser system described above. 

Since the laser system provides a measurement of the 
mean velocity parallel to the line of sight in the coherence 
volume (or, by time differentiation, its gradient: dV|j /dr) , 
and the radar detects the rras refractive index fluctuations, 
simultaneous measurements at comparable sensitivities of 
scattering from the upper atmosphere can bo expected to be 
a fruitful source of infoannation on atmopsheric motion and 
turbulence. A laser system such as the one described cannot 
be considered as an operational tool for detection of high 
altitude CAT on a routine basis because of the constraints 
of clear weather and good seeing conditions. However, as a 
research tool, it appears to offer considerable promise if 
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TABLE II 

Comparison of Siqnal/Noise per Pulse for Ground-Based Laser 


and Radar Detection at 10 Icm Altitude 


System Characteristics 


Wallops Island 

Radars 

10 . 6 Coherent Laser 

Pulse Energy 

Resolution 

Pulse energy: 200 mj 

UHP 6 
S band 0 
X band 1 

joules 

..8 

150 m 
150 m 
300 m 

H OSS'-* 

Aperture: 50 cm 

Efficiency: 0=0.02 

Absorption loss: 6 db 


Signal/Noise Ratio per Pulse 


Radar 


2.1x10 

0.21x10 


Laser 


no (10 . 6p) 


10 ^ km-^ ster ^ 


(S/N) 

(single 

pulse) 


S W 

(600 pulses- 
Incoherentl 






















an overall efficiency factor .of 2% can actually be achieved 
and if the backscattering conditions at 10. 6y are as high 
as assumed here. 

In general/ two modej of operation may be envisaged, 
one where the line of sight is fixed in space and the 
returning signal is spectrally analyzed or filtered to 
yield the distribution of velocities in the scattering 
medium, and another mode which requires measuring only the 
mean frequency offset or mean velocity of particles in the 
coherence volume and presenting this data as a function of 
range and angular inclination of the line of sight. This 
latter mode of presentation has the advantage of allowing 
presentation of the local mean velocity as a function of 
range and elevation angle while only being about 3 db in 
signal/noise belov/ the full bandwidth signal. 


..J 
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TII. MOUNTAIN LEE WAVES - A REVIEW 

Mountain lee waves are an important source of CAT. 

An authoritative and comprehensive survey of mountain lee 
waves theories/ including many historical references and 
selected references of, recent work in the USSR, can be 
found in Miles (1969) . A more recent review v/as given by 
Vergeiner (1971) . We shall draw from these and other 
sources to describe the phenomenology of the mountain lee 
waves and then discuss the implementation of a linear model 
to represent the velocity fields. In Section IV, we discuss 
the nature of a simulated LDV sensor response. 

A. Phenomenology of Mountain Lee Waves 

Standing wave trains of wavelengths of the order of 
15 km and vertical amplitudes of perhaps .5 km are frequently 
found in the lee of mountain ranges and may extend for 
several complete cycles downwind. Many v/ell-documented 
cases are given by llolmboe and Klieforth (1957) and Aanensen 
(1965) . This stationary pattern of waves may also extend 
with substantial amplitude upwards into the stratosphere, 
v/here it is occasionally made visible by mother-of-pearl 
clouds in the height range of 20 to 30 km (Hesstvedt, 1958) . 
In general, the vertical wavelengths are of the order of 
12 km and horizontal wavelengths of 50 km or more; the 
periods (the inverse of Brunt-Vaisala frequency) arc in 



. ! 



the range of 10 to 200 minutes. The propagating velocity 
relative to the wind at cloud height is about 80 m/sec. 

The ir.ouiitain lee waves are basically a regulated motion of 
a stably-stratified atmosphere when it is perturbed. The 
perturbing force pushes the atmosphere away from equilibrium 
while the gravitational force acts to pull it back. Much 
of the phenomenology of the lee waves can be visualized by 
considering a stable configuration composed of two flowing 
layers of different density. As they flow over a finite 
mound, the fluids are displaced so that part of the total 
energy is transferred to potential energy to move the heavier 
fluid upwards. However, due to the conservative nature of 
the system in the absence of dissipative mechanisms, a con- 
tinuous interchange between potential and kinetic energy is 
established so that an oscillatory motion (internal gravity 
wave) is formed. Depending xipon the nature of the fluid 
layers, their degree of stratification, the vertical dimen- 
sions of the fluid, and the I'elative dimensions of the mound, 
it is conceivable that certain modes of the oscillation 
(resonant waves) are more excited than others. In the lee 
of the mountains , these resonant waves are the only ones that 
do not decay rapidly v/ith distance downstream and are usually 
referred to as the 'lee waves*. 

We can consider the surface of discontinuity between 
these tv/o fluids mentioned above as a wave guide on which 
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the gravity v/ave travels; as a matter of fact, ocean surface 
v/aves are another extreme example of ducted waves. In the 


atmosphere, hov:ever, no such discontinuous surfaces exist, 
but layers of varying stratification (or layers of different 
stability 3 = ” ^ are imbedded within each other. The 
layer of higher stability (S large) imbedded in a region 
of lower stability (3 small) forms a duct or wave guide 
for horizontal propagation of internal gravity waves. 

This phenomenon was demonstrated by Whitham (1971) , 
starting from the dispersion relation for the gravity wave 


2 

(0 




( 1 )^ is the Brunt-Vaisala frequency. Whitham indicated that 
no gravity waves are possible except when w £ Uq. But 
because Vp x V (|-) = g ^ x $ = 3, where 0 is potential 
temperature and k is now the vertical unit vector, the 
amount of vorticity generated in a stratified flow is 
exactly the square of the Brunt-Vais^l'A frequency. There- 
fore, the layer of higher stability can carry higher 
amounts of vorticity or disturbance than its neighboring 
layers. If the vorticity contained in the gravity wave 
of frequency o) is greater than either viscous damping 

or acoustic radiation must take place in order to main- 
tain the wave form. In fact, this was the early version 
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of Richardson’s criterion of stability^ i.e. 


or 


Q2 < S. M 
“ 0 dz 


R. > 1. 

1 


The more comraoniy adopted critical Richardson's number 
is 1/4 , which can be derived by balancing the potential 
energy required to move the heavier fluid upwards and the 
kinetic energy contained in the flow motion [see Ludlam 
(1967) 1 . Whitham also indicated that the gravity wave is 
transverse, i,e., the group velocity 

2 

= f 

'1 




Ill k k„ 
o X z 


is perpendicular to the wave vector (k ,k ) , so that a 
linear model can be applicable to leirger amplitudes than 
v/hen it is applied to the problems of (say) an acoustic 
wave. 

The importance of the ducting phenomena to the 
present study is revealed in the following aspects: first, 

only ducted modes are apt to be found at or near the dis- 
turbance source altitude when observations (such as the 
LDV measurements) of internal gravity v;aves are made at 
some lateral distance av/ay from the source. Second, since 
the ducted modes are confir, ed in the duct for a long period 
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of time, large wave amplitude resonance, for example the 
lee V7aves, often are realized. Specific studies of this 
ducting process were made by Jones (1972) and Danielsen 
and Bleck (1970) . 

The relationship between the presence of the ducted 
internal gravity waves to CAT has long been sought in hopes 
that it may explain why CAT appears as sporadic patches in 
a free atmosphere isolated from any visible vertical con- 
vective activity, and how and from where CAT derives its 
energy. Work by Scorer, Bretheren and Hines has illuminated 
this relationship. The phenomena of wave resonance, the 
ducted waves process and the discovery of the existence of 
a critical level in the atmosphere and the energy trans- 
ferring mechanism between the mean wind shear and the inter- 
nal waves are all involved. Booker and Bretheren (1967a) 
found the critical level (which is defined as the altitude 
where the wave's horizontal phase velocity is equal to the 
mean wind velocity). Above v/hich the solution shows no dis- 
turbance; however, large amplitude perturbations build up 
underneath the level so that the breakdown into turbulence 
is possible. Evidence of this critical level was presented 
by Gerbier and Berenger (1961) who reported turbulence 
found by glider pilots at such heights with little or no 
lee wave action above. The effect of this level was 
analyzed through a series of papers by Bretheren (1966, 
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1967a, 1967b, 1968, 1969a, 1969b) and Hazel (1967). In 
summary, the mechanism of this critical level is to re- 
absorb the internal gravity wave energy back to the mean 
wind flow. Bekofske and Liu (1972) demonstrated that the 
mean wind shear is enhanced by this mechanism so that the 
local Richardson number is reduced to a value smaller than 
1/4. The process eventually leads to Kelvin-Helmholtz 
instabilities and thus causes turbulence production. 

Other physical aspects of lee waves are: 1) upstream 

propagation of disturbances in a stratified medium; 2) the 
effect of the earth's rotation upon the flow; 3) three- 
dimensionality effects. 

The disturbances caused by a barrier in an ordinary 
homogeneous fluid flow usually decay very rapidly upstream. 
Therefore, the description of the upstream boundary condition 
is not very critical. However, in a stratified fluid flow, 
the disturbcuxces propagate far upstream by riding on the 
internal waves , so that the details of the upstream con- 
ditions are quite important. In reality, viscosity ensures 
that the disturbance is dissipated at a finite distance 
upstream. But unless the practical application v/arrants 
the mathematical complexities required to include either 
viscous effects or to treat the problem from an unsteady 
point of viev/, it is necessary to remove the upstream 
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perturbations a priori . Scorer (1949) arbitrarily added 
terras to the solution to cancel disturbances generated at 
positions far upstream; Danielsen and Bleck (1970) and 
Vergeiner (1971) chose the complex integration contour to 
ensure that the upstream disturbances are very small. 

For a small enough wavelength the effect of the 
earth’s rotation can be neglected. This can be shown by 
comparing the quantities gli to The former represents 

the maximum potential energy required to flow over the 
mountain height H ('vl km) ; the latter measures the kinetic 
energy due to the earth's rotation. Substituting typical 
values in lower atmosphere, U (wind speed) '\>lQ0m/sec, f2'v7xlO 
rad/sec, X'v-50 km, v/e have 35. Therefore, unless 

waves have wavelengths of thousands of kilometers (i.e., 
comparable to the earth’s radius), the energy derived from 
the earth's rotation contributes little to the wave motion. 

Although a fully three-dimensional analysis can be 
formulated much the same way, for example by Scorer (1953, 
1956) , nearly all analyses in the literature are two- 
dimensional. This is partly due to the dearth of meteorolog- 
ical data collected on a horizontal plane and partly because 
mosf. essential information can be derived from the under- 
standing of a simpler two-dimensional analysis. Basically, 
Scorer's investigations indicate that in addition to the 
waves on the vertical plane containing the flow direction. 
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transverse wave systems are created by the lateral flov/ 
around the mountain. The disturbances behind the mountain 
v/ould interfere with each other. The portion of the air 
flow around the mountain reduces the mass of air passing 
over it and leads to a reduction of v/ave amplitude and a 
comparatively rapid damping downstream. Therefore, one 
can believe that a 2-D model should provide an upper bound 
upon the wave amplitude. One should bear in mind, however, 
that three-dimensional effects, such as the convective 
instability and the preferred ducting cones (Jones, 1972) 
are excluded. 

The present model will be formulated according to 
the following assumptions: the processes are rapid enough 

to be adiabatic, the disturbances of velocity are small 
compared to the free stream mean wind, and the displacement 
is a small fraction of the layer depth. The earth's rotation 
is negligible, and the Boussinesg approximation is valid 
for gravity v;aves in a compressible atmosphere. 

We will follow a procedure here similar to that of 
Danielson and Bleck (1970) . The stratified atmosphere is 
divided into several layers; each is assigned a constant 
value of Scorer’s parameter. The linearized equation of 
motion is Fourier-analyzed in the horizontal direction. 

The resonant v/ave modes are then solved as an eigenvalue 
problem. The inverse Fourier transform is performed on 
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the complex plane as a contour integral. Contours are 
chosen to include the poles and to minimize the computation 
Sample calculations are performed for two sets of data 
collected by Lilly (1971) on February 18, 1970/ along the 
cross-section betv;een the Kremmiling and Colorado VOR 
stations/ and on February IS, 1968/ near the continental 
divide (given by Lilly and Toutenhoofd (1969)). 








B. Forinulation of the Model 
1. Equation of Motion 

The linearized equations of motion for steady horizontal 
two-dimensional flow of a vertically stratified atmosphere in the 
Boussinesq approximation have been shovm (Bretherton/ 1966 
and Vergeiner^ 1971) to be reducible to the form 

2 2 

^ + K(z)4> — 0 (1) 

v/here x is the direction of the unperturbed v;ind and z the 

2 

vertical direction. The coefficient K (z) is given by 


K^{z) 


a M _L 

0 3 2 y2 


1 

u 


3^U 

3z^ 


( 2 ) 


Here U(z) is the horizontal wird and 0 is the "potential" 
temperature : 



where T is the temperature at altitude z, p(z) is the pressure 
at that altitude and p^ is a reference pressure (conveniently 
the surface pressure) . The quantity ({> can be taken as the 
streamline vertical displacement, the streamline slope or the 
vertical velocity as desired.- We will find it convenient 
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to choose (j)(z) to represent the streamline displacement since 
this quantity is continous even though there may be dis- 
continuities in the horizontal wind profile U^{z). 

Several investigations have shovm hov/ lee waves develop 
dovmstreain of mountain ridges {Scorer, 1948; Foltz, 19 66; 
Danielsen and Bleck, 1970) . We follow a procedure here similar 
(although not identical) to that of Danielsen and Bleck. We 
think of the atmosphere as being divided into a number of hori- 
zontal layers in each of which we approximate K(z) as having 
a constant value (Figure 5) . The ntli layer lies between the 
altitudes of z„ and z and in this layer K(z) = K . 

To solve Equation (1) subject to boundary conditions 
at the bottom and the top of the atmosphere, we Fourier analyze 
the lower boundary profile in the form 


h(x) 



ik X dk 

X X 


(4) 


The general solution for the displacement may be written in 
terms of its Fourier components. 


<i> (z,x) 



ik^x dk 
e 


( 5 ) 


The equations for the coefficients <j). have the form 




J : 1 




92 , 


^ + (k^(z) - = 0 


(6) 


In each layer we assume K(z) to have a constant value 
(kj^) . The solution in the nth laj r is 




+iK„(z-z„) 


-itc (z-z ) 




( 7 ) 


where 
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n X 


Continuity of value and slope* relate and to 


and 


XK 
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An e ( 


^n+1 " g g"^-'^n(^n+l “ ^n) 


•^n+l r^n+1 ®n+l 


In terms of the vectors 


, = C 4> 
^n+1 n ^n 


*n '(b") 


we have the recursion relation 

( 9 ) 


v/here c (k ) is a 2x2 matrix of the form: 

n X 


Continuity of slope is required to conserve mass 
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(10) 


2.1 The upper Boundary 

At the top layer v/e have two possible boundary con- 
ditions. For the first, the level is taken to be a 

rigid horizontal surface (i.e., there is no motion in the 
Nth layer). Here v;e must set A^j ^^N+1~^N^J 

+Bjj ^^N+l”^N^J ~ other boundary condition 

we choose the Nth layer to be infinite in extent and require that 
no energy be propagating downwards from above. This state may 
be achieved only by setting B^j = 0. Thus our upper condition 
has one of the following forms 


or Bjj = 0 • 


2.2 Lower Surface Boundary Condition 

For the general case of a homogeneous upper boundary 


condition we can write 





V7here, in our tv'o cases, 


= (J) (Lp [-2i 




Using the inverse relation to Eq. (9) 


lj> = C '1' lT t 
^n n ^n+1 


we can write 


\ ^N-1 ^o 


Thus v/e may deduce by repeated application of Eq. (13) 


'i’l = s-1 n^i 


It is convenient to define the 2x2 matrix according to 


C^_j_ C^_2 • • • • ; 


thus 


'> ■ C;) = 


A„ , e 
N ^N“l o 


The vector c„ t e can be written in terms of two parameter 
N-1 o 

C-(k„)and C„ (k„) according to 


Now our lower boundary condition says that 


Ck 


(19) 


X 


where is the fourier component of the surface displacement 

^x 

Thus 

( 20 ) 


X 


= VCg 


Denoting the diagonal function matrix as 


iK^(z-Zn^ 


P ^ 
n 


( 21 ) 


-iK„(2-zp 


we can write the solution (for the k component) in the form 


X -1 

'*'^x ” *^n~l ^N-1 ^o ^n ® 


( 22 ) 


and thus the general solution in the form 


(z,x) = 


Cl 


X -1 

C^+Cg *^n-l ^N-1 ^o n 


F„ dk^ . (23) 
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Here ^ ^ ^ ? ^nd P^ all depend on k^ . 


Resonant modes (those waves that can propagate far 


from the region of the disturbance: lee v/aves) are those 

free modes for which the integrand has poles. These will 



1 


correspond to the zeros of the factor (C^+C^) along the real 

axis (these zeros will always occur in pairs at +k ) since 

. 2 
K involves only . 

In order that the integral in Eq. (23) have a unique 

value, it is necessary to specify on which side of the various 

poles the integration path passes. This is easily determined 

by the following considerations. Suppose that we add a small 

but finite amount of viscous damping to the atmosphere. The 

poles will then be shifted off the real axis by a small 

amount y ( k ->-( k j I resonant contributions 

y ^i V i/real J 

(corresponding to the residues of Eq. (23) v/ill have an x de- 
pendence of the form 


-yx + 
e 



Viscous effects will result in positive values of y. Thus, 
taking the limit of y->0 (to the inviscid case) , we find that 
the integration path in Eq. (23) passes below all the poles 
of the integrand on the real axis. 

The direct numerical integration of Eq. (23) along 
the real axis is very difficult because of the presence of 
the poles, and it is convenient to use the theory of residues 
to obtain a more convenient path. Consider the integration 
paths shown in Figure 6. For x>0 it is necessary to choose 
a path constrained to the upper half plane in order that 
convergent integrals result. Plere we may write 






J \ J 


(residues) 


poles 


^ J 1 

2 



(24) 


On the other hand, for x<0 , a path constrained to tlie lower 
half plane is appropriate. Here 

/ [ ] = / [ ] 

1 3 


The integrals over the paths 2 or 3 represent the non-resonant 

contributions to the wave motion. In general, they arc symmetric 

with respect to the mountain and are significant only in the 

neighborhood of the disturbing boundary due to the factor 
ik X 

e ^ in Eg. (23) . On the other hand, the sum over residues 
give rise to the lee waves which can be significant well 
dov/nstrcam of the boundary disturbance. Integration of the non- 
resonant contribution of Eg. (23) can further be simplified 
by using the fact that the integrands evaluated on the contours 
2a (or 3a) and 2b (or 3b) are complex conjugates to each 
other, such that only half of the contours need to be taken, 
and the result is the real part of (2 x integrand of (23) 

X dk ) along contour 2b or 3b. The present choice of the 
contours greatly simplifies the integration of the non-resonant 
contribution without generating any artificial discontinuity 
at the mountain peak (Scorer, 1949) or requiring the cumbersome 
operations pursu( 2 d by Vergeiner (1971) to remove the singu- 


larities . 



3. Transformation to Physical Variables 

The perturbed variables u' , w* , p’ , p' / T* and 0' 
can be obtained by simple derivations from the continuity 
equation and the kinematic condition. Considering the flov7 
confined in a stream tube which is formed by the streamlines 
(j)(x,z) and c()(x,z+Az) and knov/ing that the fluid is incompres- 
sible, we can show that the horizontal velocity is given by: 


iz = (U+u') ^ Az + (j)(z+Az) - <j>(z)^ 


ii 

sz 

l+M 

3 z 


where Az is the spacing between the same streamlines at x = -to. 
By kinematic condition, we have the relation for the vertical 
velocity 


w' = U 


Siabstituting these relations to the incompressibility condition. 




we have 


1 ^ 
p dz 
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Taking into account the adiabatic condition and the energy 
equation, we have 


£' = „ 1 |£ (j, 

P P dz ^ 


and using the ideal gas law we will have 


if " I £)♦ 


1 ^ 
T 3z 


p dz/*^ 0 


1 30 A 
0 9z 


From Eqs. (26) through (31) , one can then derive the local 
Richardson number. By assuming that the flow will become 
turbulent when it is unstable (R^ < k) and that irregularities 
will continue to exist until the local Richardson nvimber 
increases to 1, we can locate regions in v/hich we expect 
CAT to exist in the lee of the mountain. The Richardson 
number describes conditions under which the flow becomes un- 
stable, and indicates where turbulence is likely, but it 
doe not allow estimates to be made of such important turbu- 
lence characteristics as the magnitude of the velocity fluctua 
tions and their length scales. This information would have 
to be extracted by other means and v;ill be disucssed else- 
where. 


These investigations are divided into two parts; a 
three-layer test case, and two seven-layer cases chosen to 
represent observations. The first part was undertaken to 
prove the convergence of the present linear analysis, to 
reveal the effects of the different boundary conditions 
applied on the top, and to test the present computer code. 

In the limit as the thickness AZ 2 of the middle 

layer approaches zero, with the product of thickness and 

2 

Scorer’s parameter K 2 AZ 2 being constant, the resonant 
wavenumber found by the present numerical scheme converges 
to that given by the relation derived in Appendix III. The 
close agreement between numerical and analytical results 
gives us confidence that this code is working properly. 

Boundary conditions for this problem are not uniquely 
defined. Solutions to Eq. (1) are presented in Figures 8 
and 9, with the input vertical profiles of the mean wind and 
the potential temperature shown in Figure 7 . Figure 8 shows 
the displacement profiles when the solid wall condition 
is applied on the top boundary; a purely real resonant mode 
corresponding to wavelength X = 118.5 km was found and the 
v;ave amplitude is about 1.5 km. Figure 9 gives the result 
when the evanescent condition is applied; one resonant 
v/avenumber was obtained with a small imaginary part, and the 
wavelength is 68.87 km. The smaller wave amplitude in tho 
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middle layer in Figure 9 indicates the effect of energy 
leaking through the upper boundary, i.e., the waves are not 
completely ducted. The continual radiation of energy up- 
V7ards results in the decreasing of wave amplitudes down- 
stream even without friction, in contrast to the case of 
a rigid upper lid where the energy is reflected back and 
remains in the floi'?. Notice that in both cases the lower 
boundary condition is satisfied. 

Neither of these tv/o upper boundary conditions is 
exact. Evidence exists that the upper atmosphere extracts 
energy from the lower troposphere (this corresponds to the 
evanescent case) . On the other hand the solid wall case 
may approximate the critical level. The existence of a 
more realistic upper boundary condition for the linear 
formulation is still in doubt. One may avoid the diffi- 
culty by extending the top boundary to infinity, but then 
the flow cannot be treated from the adiabatic and the 
linearized point of view, nor does there exist any obvious 
physical boundary upon which the above conditions apply. 
Nevertheless, the critical level formulation by Booker and 
Bretherton (1967a) represents a realistic upper boundary 
condition and this is approximated by the solid v;all 
condition. One should notice, hov/ever, as will be shown 
belov/, the solutions obtained by applying either the solid 
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wall or the evanescent upper boundary condition are in 
raany cases nearly indistinguishable. As long as the wave 
amplitude is small/ v;hich boundary condition applies is 
not a critical question. 

The second part of this investigation deals with a 
number of real cases. The vertical profiles of the mean wind 
and the potential temperatures are shown in Figure 10a and 
lla, which approximates the data obtained by Lilly (1971) on 
February 18/ 1970, along the cross section between the 
Kremmling and Colorado VOR stations, and on February 15, 

1968/ near the continental divide given by Lilly and 
Toutenhoofd (1969). In contrast to the K(z) profiles used 
in the first part of the computation, the K(z) profiles for 
these cases, as shown in Figure 10b and 11b have a minimum 
in the middle layer. The air is less stable close to the 
ground and is increasingly stable upwards such that one can 
expect that most of the CAT should occur beneath the stratos- 
phere . 

It is generally found to be true that the strong wind 
shear and stable layers are the two factors conducive to 
gravity v/ave formation. To manifest this, the chosen two 
cases for which lee wave activities were observed are espe- 
cially interesting. Actual computation for both cases was 
difficult, as it turned out, due to their extreme properties. 
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Figure 10b, Scorer's Parameter for Motcoroiogical Data on 
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In one case the Scorer parameter K (z) in the ground layer 
was large and in another it becomes negative. The first 
case occurs when the ground layer is very stable and the 
wind shear is not strong; this usually occurs in the 
evening (the February 15, 1968, case) such that the flow 
is blocked upstream of the mountain, as was observed by 
Lilly and Toutenhoofd (1969) . In this blocked area the 
perturbed horizontal velocity must be comparable to the mean 
flow in order to produce a stagnation point; in consequence, the 
small perturbation assumption breaks dovm. Another way of 
looking at this phenomenon is that blocking occurs when 
the gravity wave fronts are vertical; i.e., the disturbance 
will propagate to upstream infinity and hence contradicts 
the definition of the mean flov/ condition. Unless a priori 
information is obtained about the blocked profile, a detailed 
solution in this region cannot be expected from the linear 
theory. The second case corresponds to a nearly neutral 
stable ground layer with strong wind shear; this usually 
appears in the daytime profiles (for example, the February 18, 
1970, case) so that the flow can easily be churned into 
vortex motion. 

To avoid dealing v/ith imaginary K(z) values, which upon 
substitution into Eq. (7) will lead to exponential solutions, 
Foldvik (1962) , Conover (1964) and Lovill (1969) simply 
dropped the second shear term in Eq. (2) . Similarly, 
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Danielsen and Bleck (1970) applied heavy smoothing to the 
2 

K (z) profile, and Vergeiner (1971) also applied discon- 
tinuous shear profiles. However, this shear term affects 
the vertical convection of the horizontal mean momentum and 
the horizontal convection of vertical perturbation momentum, 
and should be included in evaluating the flow. (A physical 
argument could be advanced to remove this difficulty, but 
was not applied in the present study: It is well known 

from Rayleigh’s criterion for an inviscid unstratified flow 
that a flow with a point of inflection in its velocity pro- 
file is unstable; it can be shown that this criteria still 
holds for the stratified flow. For the wind profiles in 
Figure 10 for February 18, 1970, there are two points of 
inflection, turbulence must ex'‘st to reduce the wind 
shear and restore the flow staLility. Without dealing with 
a nonlinear problem, one could be content v/ith the assump- 
tion that the flow there should be well mixed so that no 
gradients exist in either velocity or in the potential 
temperature profiles. Therefore, one could assign a small 
positive value to K (z) wherever it is negative, and the 
solution V7ould not be found to be sensitive to the assumed 
value, as was shown in the first set of computations.) 

Difficulties were encountered due to the lack of 
information in the ground layer, represented by the dotted 
region in Figs. 10a and 11a, from 3 to 5 km above sea level, 
even with the best available meteorological data. Unfortunately 



leewave solution depends critically upon the input mean 
flow, the wind shear profile in particular. To demonstrate 
this. Figure lib shows the three possible K (z) profiles 
near the ground for the February 15, 1968, case; the real 
profile could be anyx^here between them. The results are 
shown .1 Figure 12, In cases a and b there is little wave 
action, while case c produces a leewave pattern with some 
similarity to the experimental data given by Lilly and 
Toutenhoofd (1969) - The resonant wave numbers are tabulated 
in the following table. 


resonant 

wave 

numbers 


Real 

Imaginary 


Real 

Imaginary 


,1315699 

,1132932 

.573351 

,0213953 


Case 


b 

c 

.14943 

.228069 

,081939 

.028817 

.573387 

.573475 

.021373 

. 0213169 , 


The second resonant wave number does not vary significantly 
as it is not controlled by the ground layer property. It is 
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interesting to see that as th.e ground K'(z) value gradually 
increases from negative values, the imaginary part of the 
first resonant wavenumber gradually decreases. Figure 12a 
shows rapid decay dov/nstream and increasingly visible 
wave activities in Figures 12b and 12c, 

Both cases investigated were really pathological 
cases for the linear theory. However, one call still get 
realistic results if one is willing to give up the flov/ 
details in the ground layer. Linear adiabatic theory is 
not valid near the ground due to the existence of the tur- 
bulent Ekman planetary boundary layer, blocking, and non- 
isentropic processes such as viscous damping and heat 
transfer. Here we only attempted to preserve infoannation 
far from the disturbances without invoking trial and 
error methods for fitting the blocked profile or applying 
nonlinear boundary conditions. 

2 

The Scorer’s parameter k (z) determined from the 

dotted profiles in Figure 11a near ground is also shown 

in dotted lines in Figure 11b. It is difficult to assign 

2 

an arbitrary constant value for k (z) in the ground layer, 

since it varies from 2.25 to -.25 over a height of only 

2 km. However, if we take the same case and ignore the 

small scale structures (variations over dimensions less 

2 

than 1 km) and recalculate the k (z) , wc obtain the dashed 
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. Sensitivity of F3.ov/ Solution Upon the Lower 
Portion of Meteorological Data 


Figure 12 
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Figure 12. Sensitivity of Flow Solution Upon the Lower 
Portion of Meteorological Data 






lines in Figure 13. The solid line in Figure 13 represents 
the approximate profile actually applied in the calculation. 
The results are showii in Figure 14 through Figure 16. 

Figure 14 shows the streamlj .ne pattern where lee wave 
motions are apparent; all scales are in km. Hatched areas 
represent the cross-section of topography over the contin- 
ental divide. A velocity vector plot is shown in Figure 15, 
where strong up-and-down motion over the mountains can be 
seen. Figure 16 shows the distribution of Richardson's 
number. The dots represent unstable regions (where < 1/4) 
and stars show where regions are marginally unstable. 

(1/4 < < 1) Unstable and marginally stable regions 

can be identified as those where CAT is likely to be 
found. 

In Figure 16, several horizontal lines at different 
altitudes and three lines at angles of 6.6®, 45“ and 171.5° 
are drawn to simulate the lines-of-sight for an LDV system. 
Along each line-of-sight , the parallel velocity, its 
gradient along the line-of-sight, potential temperature 
and its gradient and Richardson number can be calculated. 

A series of these plots are given in Figure 17 for 0=45° , 
where 0 is the elevation angle. From these types (, j. p^.^ts, 
Table III was compiled for different lines-of-sight. om 
these data we will establish, in Section IV, the basic 
requirements on range resolution, integration time, and 
maxiiuum detectable range for a possible LDV system. 
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Figure 13. K ( 2 ) for the Meteorological Data on 

February 15, 1968, over the Continental Divide 
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Figure 16 . Distribution of Richardson Number 
(Note difference in vertical and horizontal scales) 
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Figure 17b. Gradient of parallel velocity along line-of-sight 
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Figure 17e. Richardson number along line-of-sight 


Figure 17. 


Typical Inputs for LDV Sensor Response 
Simulation. 





a 

Altitude 

0® 

3.7 km 

0° 

5.7 km 

0° 

9.7 km 

wm 

ground 

45° 

ground 


Maximum* 
3v ii/Dr 

Range at 
Maximuia 




17 km 


4 km 


8.5 km 


8.5 km 


8.5 km 


4.25 km 


1.3 km 


Table III. Compilation of Atmospheric Wind Data 

The temporal variation of wind velocity cannot be 

obtained from the lee-wave model, but a rough estimate can 

be made here. Tv;o time scales are relevant here: a long 

time scale associated with the wave motion in a stratified 

(a 

medium { ^ and a short time scale for dissipation 


of small eddies /v/g , where V is kinematic viscosity and 
c is the dissipation rate. The former is in the order of 
several tens of minutes, while the latter is in the order 
of .01 second. Any time scales between those two extremes 
are possible and are associated with medium size eddies. 


































IV. ASSESSMENTS OF CURRENT SYSTEMS 


For a pulsed system , we have constructed Table IV 
from the results in Section III. This table relates pulse 
length to vertical range resolution, peak AVjj 

for several elevation angles . 


Elevation 

Angle 

Pulse Dur- 
ation T 

Pulse 

Length 

Ct/2 

Vertical 

Resolu- 

tion 

AV , 
pulse 

=X/2t 

Peak j t = 

(3v‘\ cT 
) * 2“ 

0 

o 

11 

6.7p sec 

1000 m 

0 

.79 m/se 

~4 m/sec 

2.2n sec 

330 m 

0 

2.4 ra/se- 

r ~1,3 m/sec 

CD 

It 

o 

6.7u sec 

1000 m 

114 m 

.79 m/sei 

3 "4 m/sec 

2.2y sec 

330 m 1 

38 m 

2.4 m/se< 

3 -“1.3 m/sec 

0 = 45° 

6.7ii sec 

1000 m 

707 m 

,79 m/se< 

3 ~27 m/sec 

2.2y sec 

330 m 

230 m 

2.4 m/se^ 

3 "“9 m/sec 


Table IV. Characteristics of a Pulsed System at Various 
Elevation Angles 


At lower angles 0= 0°, 7°, higher range resolution 
can be achieved; however, at low angles, it takes longer 
ranges to probe higher altitude wind information. At 
6= 45° , the range resolution is poor but higher wind vari- 
ation exists along the line-of-sight. 

To assess the measurability of turbulence level by 


a pulsed system, it is important to compare the root mean 
square turbulence fluctuation expected in CAT to . 





From the available data collected by Vinnichenko, et al. 
(1903), we have plotted V^^^ versus pulse time in Figure 18 


The V 


is obtained by 


c — 
2 k 


and k = ct, C = 1.375. Four levels of turbulence are given, 
each corresponding to a different dissipation rate which 
is directly related to the severity of CAT experienced by 
commercial aircraft pilots . se ” ^ also shown in 

Figure 18. It defines turbulence levels measurable for each 
pulse time. For example, the 2^1 second pulse system cannot 
measure turbulence belov7 moderate levels. 
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J. Alex Thomson and Mark F. Dorian 
1.0 INTRODUCTION 

In this appendix we present an analysis of a general 
bistatic laser doppler velocimeter system and its appli- 
cation to the measurement of atmospheric winds in the 
presence of turbulence. This analysis is excerpted from a 
draft of a report originally co-authored in 1967 by one 
of the present authors (JALT) . Although several copies 
of the report have been circulated privately, no formal 
distribution has ever been made. Part of the v;ork covered 
here has been treated subsequently by Sonnenschein (1970) — 
particularly that for the coaxial configuration in the 
absence of turbulence, and that work follows closely the 
approach described herein. However, since neither the 
bistatic configuration nor the effects of curbulence appear 
to have been treated elsewhere, we include the entire 
treatment here. Section 4.2 of this Appendix, "Atmospheric 
Attenuation", was updated by F.P. Boynton to reflect current 
atmospheric transmission models. 




2.0 ANALYSIS 


2.1 General 

In atypical optical heterodyne system for detection 
of scattered light, an optical transmitter projects a beam 
of coherent light into a particle cloud. After passing 
through the receiver optics {which are focused at some 
particular point in the incident beam) , the scattered light 
is mixed with a coherent reference beam. The mixed beam is 
then detected by a square law photodetector and the spectrum 
of beat frequencies is diagnosed to evaluate the particle 
velocity distributions. The transmitter optics may act 
also as the receiver optics (for direct backscatter measure- 
ments) or separate optical systems may be used. The antenna 
properties of such detection systems are well understood 
in the radio and microwave frequency region and, as is well 
known, the identical analysis applies to optical heterodyne 
systems; the only differences of importance to the present 

discussion are associated with the dimensions of the near 

* 

field of the antennas . For easily available optical aper- 
tures {say up to 30 cm) , the length of the near field for 
visible wavelengths can be as large as 200 kilometers, where 
as at common radar v/avelengths , for example 10 cm, the near 

The near field (Fresnel or transitional region) extends to 
the range at v?hich diffraction has doubled the width of a | 

collimated beam. Only within the near field can focusing [ 

be achieved. I 


field for antennas up to 10 meters diameter is less than 
one kilometer. Thus, in applying results from radar and 
microv;ave theory to optical heterodyne systems, it is im- 
portant to remember that theory valid only for the far 
field may be inapplicable at typical atmospheric ranges 
for some optical systems. This point is critical to the 
present application since it turns out that the greatest 
sensitivity for detection of aerosol motion at moderate 
ranges (-01 to 1.0 km) in the atmosphere is achieved at 
direct (6 = tt) back scatter- The only way of obtaining 
range discrimination with CW monostatic systems under these 
conditions is to focus on the region of interest and this 
can only be accomplished by working in the near field. In 
the far field it is necessary either to utilize pulsed sys- 
tems or to work with crossed beams (i.e. , scattering angles 
different than it) . 

In the present analysis, the Fresnel-Kirchoff formu- 
lation is used to evaluate the signal-to-noise ratio for the 
heterodyne detection of the Doppler-shifted light scattered 
from a cloud of moving particles. However, before proceeding 
to this detailed analysis, it is instructive to treat a 
simplified model of the scattering and detection process 
that' contains the basic elements of the more detailed 
formulation but not the complex formalism. 
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2.1.1 Simple Model 

There are tvro quantities v;hich must be evaluated 
to determine the signal-to-noise ratio that can be expected 
in a Doppler-scattering measurement of a particle cloud; 
the level of the signal return from an individual particle 
and the number of particles or volume of space that can 
be heterodyned efficiently at one time. If N is the 
total number of photons transmitted during the observa- 
tion time and a the particle scattering cross section (at 
the appropriate angle' y then the number of photons scattered 
by a particle located v^ithin the transmitter beam at a 

station v/here its cross sectional area has the value 
2 2 

d^ is N c;/d^ . Assuming the area of the receiver aperture 
2 

to be irR^ and that it is located at a range L, the number 

2 2 2 

of photons detected by the receiver is n N a irR^ /d^ (4 'itL ) , 
v/here n is the detector quantum efficiency. The power 
s3.gnal-to-noise ratio (i.e., the number of received photons 
which contribute significantly to the heterodyne signal) 
is 


N n^ a TiR^ /d^ (4irL^)j Vjj 


Here n^ is the local densit,y of particles 
range of interest and Vy is the volume of 
can be heterodyned effectively. 


in the velocity 
the cloud that 


i A 

i i; 

-i 

iS ,i 


- 
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The easiest configuration to visualize is that where 
the receiver optics are focused at some point in the trans- 
mitter beam. In the normal heterodyne system, the local 
oscillator (reference beam) may be considered a point 
source. A diffraction limited virtual image of this point 
source is present in the particle cloud at the focal point 
of the receiver optics. This image is a slender ellipsoid 
with a diameter (d^) of the order XL/2R^ and a length 
determined by the depth . of focus {«» 2d^L/R^) . Any light 
scattered by a particle in this volume will have wavefronts 
that are parallel, v/ithin a fraction of a wavelength, to 
the reference beam vravefronts and will thus give rise to 
strong beat frequencies. Light scattered from particles 
outside this volume will have wavefronts that are either 
inclined to reference beam fronts or have a considerably 
different curvature and will thus not result in broad 
fringes (i.e., will not chop efficiently). The more de- 
tailed calculations presented in the next section show that, 
although less important, the region outside the diffract.'' on 
ellipsoid is not completely negligible at backscatter. 
Nevertheless, for a very rough estimate of V^, we will use 
a value equal to that volume of the transmitter beam which 
is interesected by the virtual ellipsoidal image of the 
local oscillator (see Figure l-l) . 
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For scattering at angles large compared to both R^/L 
and to the angular width of the transmitter beam, this 
volume is approximately equal to 

whereas for direct backscatter, it is of the order 

The approximation of uniform transmitter beam width and 
negligible contributions from regions outside the receiver 
diffraction ellipsoid are poor for direct backscatter and 
render this latter expression useful only when » d^. 

When the transmitter is focused on the same point 
as in the receiver, d^ is equal to AL/2R^ and the power 
signal“to-noise ratio for scattering at a finite angle 
{sin0 >> R/L) is given by the expression 

(S/H) (n N 0 i) (i-l) 

where R is effectively the smaller of or R^ 

^R = + R^^ ^ . At backscatter for equal apertures 

(0 r and R^ = R_) : 


(S/N) is N a 


( 1 - 2 ) 


These expressions agree with the results of the more exact 
calculation described in the latter part of this section 
except for numerical factors of order 0.5. According to 
Eq. (1-2) , the signal-to-noise ratio at backscatter is 
independent both of range and of aperture, whereas for 
scattereing at a finite angle, it increases with aperture ■ 

3 

and decreases with range. Siegraann indicates that the value 
n N a X/4 (i.e., essentially the backscatter value) is 
an upper limit for the signal-to-noise ratio that can be 
achieved with normal heterodyne systems. Although it is 
theoretically possible to attain signal-to-noise ratios 
greater than this limit, the optical or electronic tech- 
niques required to achieve these increases are sufficiently 
complex and awkward that, they may not be practical. 

Comparison of Eqs, (I-l) and (I-2) shows that, at 
moderately long ranges where the angle subtended by the 
aperture is small, the signal-to-noise exhibits a strong 
dependence on scattering angle. For an aperture of 10 cm 
diameter and a range of 100 meters, the direct backscatter 
signal is three or four orders of magnitude greater than 
that for right-angle scatter. These high levels can be 
achieved, however, only for scattering angles very close to 
180** (V7ithin 0.3 milliradlans for this example). 
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At direct backscatter the range resolution for a 

2 2 

CW system is determined by the depth of focus (6L *= X L /R ) 
and, in the visible for a 10 cm aperture, is of the order 
of 2% at 100 meters range and 20% at a kilometer. At 
ranges comparable to or greater than the length of the 
near field (L ~ R /X) , significant resolution can only be 
obtained by pulsing the laser or by working at finite 
scattering angles. In this latter method, the range reso- 
lution is of the order X/2R sin0 . When the scattering angle 
is always adjusted to yield a given range resolution, the 
signal-to~noise ratio exhibits the expected inverse square 
dependence on range: 

(S/N) « (n N n^a^) SL/L^ for L » . (1-3) 

In the following section v;e carry out a more rigorous 
analysis of the scattering-heterodyne measurement for more 
general conditions than the focused configuration just con- 
sidered. The calculation is essentially a more exact for- 
mulation of the above model. In Section 3 we repeat the 
calculations including the distorting effects of atmospheric 
turbulence betv/een the scattering point and the transmitter- 
receiver system. 
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2.2 Fresnel-Kirchof 1: Formulation 


In order to obtain a formulation that v;ill permit 
fairly general optical configurations to be considered 
but also allow explicit and simple evaluation of the various 
integrals that arise, we will introduce a number of simpli- 
fying ass;ainptions. First, we assiarae that the angle between 
the axis of the transmitting lens or mirror and the receiver 
axis is sufficiently small that various cosine terms can 
be set equal to unity. This will limit the scattering 
angles which can be treated to be within 5 or 10 degrees of 
direct forward or backscatter. The extension to larger 
angles will be made subsequently in a less rigorous fashion 
on the basis of the physical picture that evolves from the 
more detailed analysis. Second, we assume the radial inten- 
sity distributions in the transmitted and reference beams 
f 2 2 2 1 

are Gaussian li.e. , |i()| ~ exp{-2r /R )j and the field 

stops in the transmitter and receiver optics do not appre- 
ciably vignette these beams. This is a fairly good approxi- 
mation for well adjusted single mode lasers (see Figure i-2) 
This approximation permits explicit and simple evaluation 
of the various integrals that arise in the Fresnel approxi- 
mation. We also assume that separate beam expanders are 
used in transmitter and reference beams so that the widths 
of the two beams may be controlled independently (see 
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Figure 1-3) . Third, we neglect any depolarization effects 
and use a scalar formulation of the wave equation. 

The optical configuration considered is shown in 
Figure 1—3. At a range L from the transmitter, the ampli- 
tude of the outgoing beam at a distance r = (x,y) from the 
optical axis is proportional to 


exp i (kL-mt) 

i(j = 1/2 A. 

^ /ir L 


CQ CO 



— CO —CO 


iTrr"^ 

Xf,. 


0 


't 
(1-4) 


dx"dy 


where r" is the radius vector to a point x" , y" in the plane 
of the transmitting lens and A^ is the total transmxtted 
flux (in photons/sec) . The amplitude has been normalized 
so that 



Equation (1-4) may be shown to be equivalent to the expression 


R^A 


if, = 


XL 


2 ^ exp|*i (kL-wt + (jj^) - ( [l + 3^^ (1-?^.)^]} 


(1-5) 


where 



2 

n = 


/XL 

L/f^ 

1-i 
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and 


1 


'f'4- = 


irr 

XL 


l-Bt 


I 


At the beam splitter or mixer in the receiver, the 
amplitude o£ the wave scattered by a particle in the plane 
L but displaced an amount r from the transmitter axis 
is given by 


^ (r ' 
^sc- 


(r ' ) = i|j exp 


[‘ 4 ^ (J-4?) 


2 2tt /“>"f 

- X 


-* ^ (r-r')^ - (kL + Awt) 


(1-6) 


where Aw(=wv/c) is the Doppler shift due to the component 
of velocity of the particle parallel to the line bisecting 
the optical axes of the transmitting and receiving lens and 
0 is the backscattering cross section of the particle. Here 
f^ is the focal length of the receiver lens, n^ is a unit 
vector parallel to the receiver optic axis (which is in- 
clined at an angle 0 to the transmitter axis) and Ar is 
the separation of the tv/o lenses. Low frequency chopping 
of the transmitted beam, although desirable in practice, is 
irrelevant to '^b,e present analysis and is ignored, although 
if present it vjould reduce the average signal pov;er by a 
factor of 2. The reference beam vyhich is added to the 
scattered signal at the mixer been chosen to be a plane 


wave: 




•f i' 
"I ? 

to*-’ 


1 1 — - -I •■ 


t 


ref 


= otA^ exp^“i(i)t + i4> ” (r'~Ar)^/R^ j (I“7) 


where a is some constant amplitude factor and ^ a constant 
phase difference. We may set ({) = 0 with no loss of 


generality. The photocurrent per unit effective area induced 


2 2 

in the detector is equal to riltffgj, + '^ref ^ electrons/cm -sec 


v/here n is the quantum efficiency (electrons/photon) . In 


a well designed system, the reference signal level is to 


2 2 

be chosen so that >> • Thus, the total 


signal current is given by the expression 


^s = 


J~w(r’) n 


^'f'sc'f'ref '^sc'^ref^ electrons/sec (1-8) 


and the reference current is 


iR -lw(r*) electrons/sec 


(1-9) 


Here W(r’) is the apodization or transmission function for 


the receiver aperture and, in the present analysis, is assumed 

2 


equal to unity over the region where | 1 0 . The 

signal current per unit effective detector area is equal 
to the real part of the expression 


2 /2a 


a n A^ — 2 *^t 


AI, 


i (Amt + ^) - 




(r’-Ar)^ 


1+3^ 


R 


/{ 

id 

\3 
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where 


(i-i 


IT ^,,2 TT *"^\2 , 2 ir -)- 

^ (r-r') -_ (r'-Ar) +x"r-<=^ 


After carrying out the integration over the receiver aper- 
ture, we find the signal current due to the scattering by 
an individual particle to be given by the expression 


j t/ 2 o R. , 

(i ) = naA^ 

^particle ^ XL 


[e3,Cr-Ar')/R^]- 








X cos ( Awt + (j) ' ) 


(I-IO) 


Here and 3^ are defined in a fashion analogous to that 
for and 3|.- The term ({)' is a phase factor which depends 
on the system parcuneters as well as the coordinates of the 
scattering particle, and the vector Ar* is the separation of 
the two optic axes at the range L. If the axes intersect 
at some range L^, Ar' has the magnitude Ar - L sin9 = (L^-L) 
sinO and is parallel to Ar. 

In the near field of both the transmitter and receiver 
(3t/ » 1) , Eq- (I-IO) indicates that the signal current 

is large only for particles that simultaneously fall v/ithin 
the diffraction ellipsoid that corresponds to the image 


of the receiver local oscillator {also a point source) . 


In the far field these two images become the usual conical 
lobes rather than ellipsoids. When a single particle is 
being tracked (i.e., it is at the focus of both the trans- 
mitter and the receiver: . r = Ar' = 0, ?j- ~ ~ 

2 2 

total pov/er signal-to-noise ratio, i^ /i^ , can be shown 
to be given by the relation 

(S/N) = [n N(irR^^) [2 cos^{Ao)t + (j)‘)^ ] (l-H) 

2 

where N is the total number of photons transmitted (=A^ t) 

and v?e have assumed the noise to be only shot noise in the 

* # 
reference beam . Since the observation time t is typically 

much greater than the beat period 2 tt/Aw/ this expression 

may be rewritten in the more familiar form 

(S/N) = nN G.G . (1-12) 

^ ^ (4tiL^) 

Here G^ and G^^ are gains of the transmitter and receiver 
antennae relative to omnidirectional antennas : 

G = 8it^ ^3 
t t 

2 2 2 
Gr = /\^ 

T- 2 

The shot noise power (i ) in a bandwidth 1/x is equal to 

2 2^ 2 

2ij^/x, i.e., to n ct 7rR /t . 
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2 ^ 
and X /4fr is the effective aperture of an isotropxc receiver . 

When all the scatterers lie within both diffraction 
lobes {far field) or diffraction ellipsoids (near field) 
of the transmitter and the receiver, the signal-to-noise 
ratio exhibits the usual proportionality to both antenna 
gains and to the inverse fourth power of the range. However, 
when the particles are distributed over larger volumes, only 
those V7ithin overlap of the two diffraction lobes or ellip- 
soids can be detected efficiently in a heterodyne system. 

In this case the functional dependence on range, aperture, 
wavelength, and scattering angle can be quite different. 

Since the signal currents due to waves scattered by 
different particles add incoherently, we obtain the total 
signal power by summing the pov/er due to each particle. 

If n^CD is the number density of particles at the range L 
that give rise to a Doppler beat frequency within the band 
pass Af of the electronic filters, the total signal power 
from all particles is given by the expression 





d vol 


(1-13) 


Assuming that the particle density is uniform over the beam 
cross section, we may express this power in the form 



These gain expressions are equivalent to the usual expressions 
for radar antennas if the effective antenna areas are taken 
to be 2 ttR^ and . 
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.2^2 


2 2^4. ^2. “ on^(L)R^R^ exp 

2 n ct a^(ttR^) 




/ 


i ( 4^)1 




(r2 + ^l)l? 


dL (1-14) 


where we have defined 


‘'t = + 3t 

and 

«r = «r/[l + ' 

In a well designed system/ the noise should be only 
shot noise from the reference beam. Thus, the noise power 
in a band width 1/t is given by 

ij^ = n ttR^ /t (electron/sec) ^ . (1-15) 


In order to obtain an easily interpretable expression for 
the power signal-to-noise ratio, it is desirable to define 
a number of characteristic parameters for the optical 
system. It will be found more convenient to represent 
the results in terms of two mean focal lengths fj^ and 
where 
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The ratio is close to unity under all conditions. 

When one aperture greatly exceeds the other, and both 

approximate the value of the focal length of the larger 

aperture system. When the two apertures are equal, both 

f^ and approximate the value of the smaller of f^ or f^_. 

Also, it is convenient to define a characteristic 
* 

range L according to the relation 


L* = fy[l + (\f2/TT\R^)^] , 

depth of field parameter AL by the expression 


Xf ,/ttR. R 


ff2 _ f2 

(E?+E?) Vt *r. 


and an angular spreading parameter y by the relation 



where 

Von. = >/ 

Sdiff = W Rfc + ®r /''®A 


In terms of these parameters and assuming that the 
optic axes intersect at a range L^, the power signal-to- 
noise ratio may be expressed in the form 



exp^-|2sin^0 (L-L^) [(L-L*) ^ 
[(L-L*) 


dL 


2 + Al2 ] 


( 1 - 16 ) 


where N is the total number of transmitted photons in the 
observation time x. 

We first consider the case of direct backscatter. 
For a uniform distribution of particles, the power signal- 
to~noise . ratio in this case is given by the expression 


(S/N) = 1 M 



0 nj 






( 1 - 17 ) 


with the major contribution coming from an interval of 
length AL. In order that a significant range resolution be 
obtained, the depth of field AL must be small compared to 


the effective range L , i.e., the apertures must be chosen 
large enough that the particles are in the near field 
{\f 2 /TTR^R^ << 1). In this case the signal-to-noise ratio 
may be expressed in the form 


f2/fi| 


(S/N) = n N nj o 4 ^^^^2 , 



1 4- 

r 2 2 



2 


X T 

Lxfj^(RtR2, 

-P 


J 


(1-18) 


In Figure 1-4 we have plotted the signal-to-noise ratio as 
a function of receiver aperture for the case f^ = f^. 
According to Eq. (1-18) , the maximum signal-to-noise and 
the sharpest range resolution is obtained only when the dif- 
fraction patterns of the transmitter and receiver match, 
i.e., when both = R^(= R) and f^ = Lq) • The power 

signal-to-noise ratio then has the value 


ir 

(S/N) = n N Uf a 4 


(1-19) 


and the range resolution is equal to the depth of focus for 
the diffraction limited image of a point source: 





For scatter xng at angles which are considerably 
greater than the beam spreading angle y, the transmitter 
and receiver antenna patterns only overlap in the neigh- 
borhood of the range point L^. In this case the signal- 
to-noise ratio is given by the expression 


(S/N) 


TT 12 

= nNn^o4 v-ir 


X/if y sin6 

* 2 2 
(L^-L ) + AL 


and the range resolution 


6L - /if y J (L -L*)^ + (AL)V»^ sine 


In this finite scattering angle case, the maximum signal- 
to-noise ratio will be obtained v/hen both transmitter and 
receiver are focused on the same point: 


L ^ sine 


Here CS/N) is the maximum value for direct backseat ter 
l^n N. n^ cr For equal apertures the range resolution 

and the signal-to-noise ratio are given by the relations 


^ - X//tF R sin0 




(S/NJ = (S/N) R//if L sine 


= (S/N)* (rVxL) (5L/L) . (1-24) 

Thus, when perfect, matched optics are focused on the range 
point and the scattering angle is adjusted so as to yield 
the desired range resolution 6L, the S/N for a heterodyne 
detection system is identical to that for incoherent 
detection N n a (■frR^/4'frL^) 6L j . This will not, in 
general, be true for other optical configurations. For 
example, if collimated tremsmitter and receiver beams are 
used (fj. = f^ = “) , the signal-to-noise ratio is given by 

(S/N) = (S/N)* R//? sin0 *>/l^ + (ttrVx)^ - (l~25) 

In contrast to the focused configuration, the signal-to-noise 
ratio as a function of aperture for collimated beams has a 
maximum value when the aperture radiiis equals /XL/ir and, at 
larger apertures, falls vrall below the value for a focused 
configuration. In general, any attempt to broaden the beam 
to accept larger solid angles (this is equivalent to 
setting f^, f^ to infinity or negative) necessarily results 
in a reduction of the signal-to-noise ratio. For the 
crossed beam system just considered, the heterodyne 


signal-to-noise ratio in the near field is reduced below 
the incoherent value by the factor [l + if 

collimated rather than focused beams are used; this would 
entail a loss of about 15 db in the visible for R = 1 cm 
and L = 100 meters. 

In Figure 1-5 vje have plotted, for various ranges, 
the maximum signal-to-noise ratio that can be attained at 
a given scattering angel with a 10 cm diameter aperture. 

In general, the levels shown can be reached only with 
matched and focused optics. The range resolutions that 
can be obtained in the visible with this aperture are indi- 
cated in the figure. It is apparent from this figure that, 

_3 

for moderate range resolution requirements (6 l/L >10 ) , 

the maximum signal-to-noise is obtained very close to direct 
back (or forward) scatter and that orders of magnitude loss 
in sensitivity may result if the scattering angle is chosen 
improperly. 


(S/M)o/CS/M) IN DB 
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SCATTERING ANGLE (DEGREES) 

Figure 1-5. Maximum Attainable Signal-to-Moise Ratio as 
Function of. Scattering Angle 



3.0 ATMOSPHERIC PROPAGATION EFFECTS 


I © 


3.x Turbulence Effects 

Atmospheric turbulence between the scatter;.ng point 
and the transmitter and receiver will degrade the signal at 
long ranges. Inhomogeneities in the refractive index in 
the propagation path cause amplitude modulation ^ spreading, 
and wandering of the outgoing beam. The effect of turbu- 
lence on the returning beam may be thought of as a spreading 
and wandering of the virtual image of the local oscillator 
coupled with a fluctuating modulation of the receiver 
antenna gain. 

The propagation of coherent beams through the atmos- 
phere has been studied by a number of investigators. 
Spreading and wandering of a coherent laser beam have been 
measured by Hinchman and Buck^ over 9- and 90-mile paths 
under clear night conditions of low humidity and stable 
winds of 5 to 7 mph. The (theoretical) beam width of 
1.25 seconds of arc had spread to 8.7 seconds at 9 miles 
and 13 seconds at 90 miles. Further, the beam wandered 

7 

about over a region several times the beam width. Straub 
has measured beam excursions up to 6 neters in a 3.5 km 

g 

pathlength. Fried has shown theoretically that, for 
transmission of a modulated coherent beam from one point 
to another, little improvement in signal-to-noise can be 
obtained by increasing the collector aperture beyond a 


G 
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certain critical diameter, which is characteristic of the 
phase structure function for the pathlength involved. 

9 

Goldstein, Miles, and Chabot have measured the effective 
aperture diameter for pathlengths of 4.0 and 23.8 3cm and 
obtained results in essential agreement with theoretical 
predictions . 

An additional effect associated with random phase 
variations is that of frequency spreading of the signal. 
Hodara^^ has calculated the effect: 


<Af^> 



<An^> 


(1-26) 


v/here L is the pathlength, L a coherence length associated 

with the outer scale size of turbulence which has been 

empirically determined to be of the order - ^ (h the 

1 + 10 '^h 

altitude in meters) , An the relative change in index of 
refraction associated with the turbulent cell, also empiri- 
cally determined as approximately 


>.A 2^ nn-12 

<An > = 10 


h 

1600 

e 


/ 


and t is the coherence time, 
c 

For a wavelength X = 0.6ii, L = 10 3cra, h = 1 !cm, 

2 —12 

<An > = 10 and t^ = 1 second, the rms frequency spread 
is about a lei lo hertz. For the present Doppler scattering 
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problem, this effect may be safely ignored. 

In general, the measured time-averaged propagation 

effects are consistent with existent propagation theory 

based on the Kolmogoroff theory of isotropic turbulence 

11 

(such as that formulated by Tatar ski ) . Most of the 
existent analyses, however, apply to the propagation of a 
collimated beam from a transmitter to a remotely located 
receiver and need to be modified somevjhat for application 
to the present scattering problem. In this section v/e 
will evaluate the signal-to-noise ratio for detection of 
Doppler scattering in a manner similar to that given in the 
preceding section but will include the atmospheric propa- 
gation effects. It will be shown that the principal effects 
of atmospheric turbulence in the optical path are a reduc- 
tion in the signal-to-noise ratio and a loss in range reso- 
lution. The analysis closely follows that of the previous 
section and only the differences necessitated by the inclusion 
of propagation effects will be indicated explicitly here. 

When atmospheric propagation effects are important, 
the integrand in Eg. (1-4) should be multiplied by a factor 
exp^xt^/^")] account for the fluctuations in the amplitude 
of the outgoing beam due to inhomogeneities in the refractive 
index along the propagation path. Here is the 

-y 

perturbation to the log amplitude between the point r" in 
the transmitter aperture and the range point r due to the 






In evaluating the expected value of the signal 
2 

power (<ig >) # we need to calculate the value of the 
expression 



,r) 






where the brackets, denote an ensemble average over 

the various possible distributions of atmospheric inhomo- 
geneities. An approximate representation of this expres- 
sion can be obtained if the variations in the log amplitude 
for the outgoing and the returning paths are assumed inde- 
pendent. This approximation v/ill underestimate the turbu- 
lence effects somewhat. With this assumption the ensemble 
average may be shown to be equal to the expression 


exp ^-^5 [0(|r’-r"* j ) + | r"-r "’ | ) j J 

where (B(r) is the wave structure function. 

For the Kolraogoroff theory of isotropic turbulence, 
the wave structure functio:i has the form 


^(r) = {r/r 


( 1 - 30 ) 


where r is a parameter depndent on the range, wavelength. 
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and level of turbulence . For a plane wave propagating 
through a turbulent medium of length %t is given by the 
relation 

= 0.0581 , (1-31) 

where is the refractive index structure constant. Since 
we are concerned here with the structure function for a wave 
originating at the range point (f) rather than for a plane 
wave, the effective length z will be somewhat less than the 
actual range L. Setting jj, = L will tend to overestimate 
the turbulence effects. 

In order to obtain an explicit evaluation of the 
expected value of the signal-to~noise ratio, we will intro- 
duce one further approximation. Instead of the 5/3 
dependence for the wave structure function in Eq. (1-30) , 
we approximate ^(r) by the simpler form (r/r^)^. This 
dependence on r would result for linear structure function 
for the temperature fluctuations and corresponds to some- 
what greater correlation for small and less for large sepa- 
rations than the more accurate two- thirds law (see Figure 
1-6) . With these approximations, the integration in 

_ ... ■ 

The parameter r^ is equal to 0.314 times Fried’s parameter 
r^. We have chosen to use this parameter rather than Tq 
in order to simplify the algebraic form of the results. 






Eq. (1-29) may be carried out directly. The result for 
the expected value of the power signal-to-noise ratio 
for a single particle is 




(1-32) 


where 


IT R, 


= 


1 + 


exp 




1 + + e^(l-5^)^ 


and 


= (E^/r^) 


The factor is defined in an analogous fashion. When 
both transmitter and receiver are focused on the particle, 
the signal-to-noise is again given by the standard radar 

radii 

and must be replaced by the effective values 

^3 + ^3 and R r /yl r ^ + R ^ . Thus , as expected 

ta'^a t ra'va r 

8 9 

from the results of Fried and of Goldstein, et al. , the 
antennas gains saturate at aperture radii greater than r^. 
Values of r^ are moderately small at the visible wavelengths 
unless the seeing conditions are exceptionally good. Typical 
values of r (= 3.18 r ), measured by Goldstein, et al. , 
are shown in Figure 1-7. In Figures 1-8 and 1-9 we have 


equation ^Eq. (1-12) j except that the aperture 


l06 






reproduced Fried's calculations of the dependence of on 
wavelength, altitude, range, and turbulence strength. 

Returning to the evaluation of S/N, we find that, 
for a distribution of particles with density that gives 
rise to Doppler shifts in the band width Af, the power 
signal~to-noise ratio is given by an expression identical 
to that for no atmospheric perturbations: 






a expf-/2 (Ii-L^) ^sin ) ^+AL^ jll 

LJ: r^dL 


* 2 2 
(L“L ) + (AD 


(1-33) 


except that the values of the various parameters now depend 
explicitly on the value of r^; 


Y = -v/ 




108 


roCMETERS) NIGHTTIME 


r„ (METERS) DAYTIME 


CO 
(T| 


t. ^ i 31 k 
: ■ c, hr : 

tu* : ftl i-9i - 




tPhn-lfe- 




■o-rra- — 


uP 'iqI 

10*^ 


RANGE R (METERS) 


Figure 1-9. Dependence of r on wavelength, range, and 

° 2 

strength of turbulence for horizontal 
paths. (e.g., r^ = 10 cm for X = .63 y, 

R = 10 km, = 10 m (from Fried^) 
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AL = f. 






R^R 
t r 

2 2 
R +R 
r r 




(1-34) 


At direct backscatter the signal-to-noise ratio is 
given by (for R^ = R^, ^r " ^t^ 


where 


(S/W) = (S/N) ‘/Af 1 + (R/r^) 


(S/W) = n N n X/4. 


(1-35) 


At large apertures (R >> r^^) , the fractional range resolution 

•k 

is equal to AL/Rr and the signal-to-noise is (S/N) (r /R) . 

a. G 

when the aperture radius is adjusted to get a given range 
resolution (6L ) , the dependence of the signal-to-noise on 
L is given by the relation 


(S/W) = (S/W)*(r^^ /AL) (6L/D 

a 


(1-36) 


At large ranges the apertures required to obtain 


shairp range resolution become unwieldy. For these conditions 




it is necessary to separate the transmitter and receiver 
to obtain adequate range resolution. In the near field 
(L << TrRr_/21) , the optimum signal-to-noise and range 
resolution are obtained by focusing on the region of interest. 
For such a focused system, the signal-to-noise ratio may 
be expressed in a form identical to Eq. (1-36) : 

2 

{S/m = {S/m* ^ . ( 1 - 37 ) 

The scattering angle required to obtain this spatial dis- 
crimination is given by 


sin0 = X//n r^ (6L/L) 

a 


(1-38) 


In the far field (L » irRr_/2X) , both L and AL 
are small compared to L and the signal-to-noise ratio is 
identical to that for collimated beams: 


(S/N) = TT ^ (S/N)* 


(r//XL) (5L/L) 

a 



and 


C£ {S/m* (r//XL)(5L/L) 

a 


~ (X//? r^) (6L/L) 

d 


(1-39) 


(1-40) 
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These results may be summarized in the following 

manner. With high quality optics a signal-to-noise ratio 

in an optical heterodyne detection system equal to the in- 

2 2 

coherent value (t) N n^ a ttR 6L/4ttL ) can always be obtained , 

in the absence of propagation losses due to atmospheric 

turbulence, if the transmitter and receiver have identical 

antenna patterns {i.e., equal apertures and focal lengths) 

and if both optics are focused on the same point in space. 

In order to achieve this limiting S/N, the transmitter- 

receiver separation must be adjusted for each range in such 

a way as to provide as much overlap as possible consistent 

with the desired range resolution 6L. At a given range 

there is generally a limit to the maximijm achievable value 

2 2 

of 6L, equal to XL /R . The achievable S/N is consequently 
never greater than n N n^ c X/4. This effect is only of 
consequence for antenna apertures large enough that the 
range point is within the near field of the antenna 
(L < R /X) . For range points well inside the near field, 
the optimum scattering angle is 180^ and the range reso- 
lution is determined primarily by the depth of focus. 

Atmospheric turbulence acts to limit the effective 
useful aperture radius at long ranges to a value dependent 
on the wavelength, range, and turbulence level 
(.0581 , At short ranges, there is also 

limit to the effective useful aperture (equal to XL /<SL) 
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which is determined by the desired range resolution. Thus, 
use of aperture radii larger than the smaller of these 
values will not appreciably enhance the signal information. 
These limiting apertures are plotted in Figure I-IO and 
I-ll as a function of range, wavelength, and atmospheric 
turbulence levels. The maKimum signal-to-noise ratio that 
can be achieved with an optimally designed CW system, 
while maintaining a range resolution capability of 10%, 
is shown in Figure 1-12. 


113 




APERT 



[; Pigure I-ll. D0pGndence of Effective Aperture on Range and Turbulence Level 

[ X “ 1 0 ]j . 








3.2 Characteristic of Natural .Aerosol 


In Figure 1-13, we show typical distributions of 
aerosol sizes in continental air. The signal-to-noise 


I for the heterodyne signal is proportional to the mean value 
'^of the product nci: . 



na 


o 


(1-41) 


where is the total scattering cross section and Y^go the 

i 

backscattering efficiency. 

The backs cattering efficiency for particles is defined 
as 4ir times the ratio of the scattered intensity per unit 
solid angle at backscatter to the total scattered intensity, 
Diermendjian has calculated the angular dependence of 
scattering for haze at a wavelength of 0.7vi (see Figure 1-14). 
A backscattering efficiency of 0.15 has been evaluated from 
these data by nvimerical integration. For rough estimates VJe 
will estimate 0 by the approximate expression 


0 ^ « 2it a*^ for a > X/2Tf 


sa 0 for a < X/2 it 


(1-42) 
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Figure 1-13 


Aerosol size 
distribution 






no = 7.89x10”^^ Y/A cm' 

At 6000 A, no is 0.132 Y.per lam and at lO.Sii is 0.00744 Y per 
Jem. These values of the total scattering coefficient and 
the inverse wavelength dependence are in good agreement with 
the detailed calculations by McClatchey, et al. (1971) , for 
his clear air model (see Figtare 1-15.), For other atmospheric 

conditions, no may be deduced from the value of the visual . 

^ 

range (visibility) which is defined as 3.9/nOg,. Thus, 



O 


at backscatter 
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Figure 1-15. Attenuation Coefficients (per km) for 
Aerosol Transmittance (Absorption and Total Extinction) 
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At 10. 6y, the model clear air backscatter cross section 
is (using y = 0,15 and into 4ir steradians) 


= 1.12x10 ^ cm ^ 

10. 6y 


(1-46) 


= 0.0011 km 
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3 . 3 Atmospheric Attenuation . 


At a wavelength of 0.5y, the predominant contri- 
bution to atmospheric attenuation will be aerosol scattering 
At 10.6ii, the characteristic wavelength of a CO 2 laser » 
absorption by atmospheric CO 2 and H 2 O will give rise to 
losses somewhat greater than scattering losses under most 
conditions . 

These losses depend upon atmospheric temperature and 
humidity, and upon altitude for given atmospheric conditions 
CO 2 absorption at the center of the P20 line can be calcu- 
lated from standard formulas and line intensity data. 
Specific expressions are given below; 

k = S/tt [ 0 -I 6 km] (1-47) 

k = (S/cq) (£n 2/ir)^ exp(a^) erf (u) 

[l6-65 km] (1-48) 

k= (S/Qp) (£n 2 /tt)^ [>65 km] . (1-49) 

Here k is the absorption coefficeint, S the line strength, 
the collision half-width, the Doppler half -width 

Cq = (2kT in 2/mC^)^ (1-50) 
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J 


and a is given as 






The collision half-width is a function of both temperature 
and pressure. 


“C “ ' V^o> < V’^> ” 


(1-52) 


where for the CO 2 P20 line si an effective broadening 
pressure related to the partial pressure of CO 2 and the total 


pressure as 


P_ = 0.3 P_ + P 


. 13 15 

and n is 0.58 . From the AFCRL line compilation , 

a® = 0.072. 

The line strength is a function of temperature according 


S{T) 


Qv(T ) Q_(T ) 
Q^(T) Q^(T) 


[1 - exp(-hcw^^/k T)j 
1^1 - exp (-hciu^/k T^)j 


fhc ] 

lir ® -tt7“ J 


(1-52) 


0 
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Here and are the vibrational and rotational partition 
functions, which for CO 2 are approximately 


Q^(T) = [1 - exp{-960/T)] [1 - exp(-1921/T] [l - exp (-3380/T)] 

(1-53) 

and 


Q^(T) = T/0.561 


(1-54) 


The AFCRL compilation^^ gives S(296°K) as 5-73xlo“‘® 

(cm^-atm) ^ ^ and E" as 1552.046 cm while the central 
296°K 

wavenumber of the P20 line is 944.19 cm 

H 2 O produces spectrally continuous absorption in the 

8-14 p region. This absorption is produced by the far wings 

of strong H 2 O lines in the 6.3y or rotational bands or 

(possibly) with H 2 O dimer formation. The absorption coef- 

o X5 

ficient for the P20 line at 296 K can be expressed as 


k = 2.50x10 ^ (P + 193 P„ ,.) cm ^ atm ^ (1-55) 

H2O 

where P and Pjj ^ are the total pressure and the H 2 O partial 
pressure, respectively, in atmospheres- No measurements 
are available at lower temperatures. The trend at higher 
temperature is that k decreases V7ith increasing tempera- 
ture. In the absence of much data on temperature dependence, 
the above expression for k can be used at atmospheric tem- 
peratures. 
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Conditions in the atmosphere vary considerably v;ith 
time and position. The U.S. Standard Atmosphere Supplements 
(1966) contain several model atmospheres v/hich can be used 
for design studies. Table I-l gives the temperature and 
relative humidity at altitudes up to 10 km for several of 
these model atmospheres. Table 1-2 gives the temperature 
structiare at higher altitudes, where the water absorption 
is negligible with respect to CO 2 absorption. The pressure 
at a given altitude can be found from the more detailed 
tables given in the 1966 tables or by using the relation 

h 

P (h) = P (o) exp J dii'/Hp] (1-56) 

o 

where 

= i ~ = .0293 T km (1-57) 

P Mg 

is the atmospheric scale height. 

13 

Yin and long have calculated the CO 2 absorption 
coefficient as a function of altitude up to 65 km for the 
P20 line center for the January and July 30®W model atmo- 
spheres. Their results for the optical depth per km of path 
length (t/H = k P„ ) are given in Table 1-3 in the form 

of polynomial fits in different altitude regions. 

15 

McCoy, Reusch, and Long have calculated the H 2 O 
absorption coefficient as a function of altitude up to 10 km 
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Table I-l; Temperature and Relative Humidity as Functions 
of Altitude for Several Model Atmospheres 


Alti- 

tude 

(km) 

15°N 
Annual 
T H 

3C 

Jan 

T H 

July 
T H 

45® 

Jan 

T H 

M 

July 
T H 

60® 

Jan 

T H 

July 
T H 

75® 

Jan 

T H j 

July 
T H 

0 

299.65 75 

287.15 

80 

301.15 

80 

272.15 

77 

294.15 

75 

257.15 

80 

287.15 

75 

249.15 

80 

278.15 

85 

1 

293.65 75 

284.15 

70 

293.65 

65 

268.65 

70 

289.65 

65 

259.15 

70 

281.75 

70 

252.15 

65 

275.55 

75 

■Bl 

— 


— 

— 

— 

— 

B 


Bl 

— 

— 

— 

Bl 

253.65 

60 

274.25 


2 
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Table 1-2: Temperature as a Function of Altitude 

for Several Model Atmospheres 

60®N 


299.65 287.15 301.15 | 272.15 1 294.15 257.15 287.15 


281.15 


286.15 

236.95 


193.15 


216.15 


203.15 

203.15 


215.15 213.15 


293.65 


203.15 

203.15 


259.15 


285.15 


261.65 



261.15 


251.15 


217.15 


219.65 


214.15 215.15 



211.15 


215.15 227.65 

233.15 236.15 219.15 238.15 


260.15 


225.15 


225.15 


249.15 

253.65 


215.15 


213.65 


207.65 


207.65 


238.65 


278.15 


271.65 


226.15 


230.15 



230.15 

237.95 


220.15 


270.15 269.15 272.15 265.65 275.65 


254.15 




269 

.15 

269 

.15 

253 

.15 

“ 

- 


265.65 275.65 


250.65 


260.15 


271.65 

277.15 


277.15 

260.15 

251.15 265.15 


241.65 


Temperature varies linearly between tabulated altitudes 

























Table 1-3: Polynomial Coefficients for CO 2 Optical Depth per km of Path 

Length as a Function of Altitude h in km (30 after Yin and Long) 


n ^ 7 ^ A 5 

Month Range h h h h h h 


January 

0-2 

6.177301“^ 

-4.243258“^ 

3.331791”^ 





2-12 

7.1772259"^ 

-1.018966““ 

5.438857"^ 

-1.440797“^ 

2.111184“^ 



12-17 

2.168201“^ 

-1.708914"^ 

4.794412“^ 

-4.821832“^ 




18-47 

2.749415"^ 

-6. 043183'^ 

4.952857"^ 

-1.73234"^ 

2.748747“^ 

-1.644807"® 


47-51 

-5.0246274"^ 

3.209868"^ 

-6.729368"^ 

4.656468“^ 



1 

51-65 

1.241526“^ 

-5.636995"^ 

8.609204“^ 

-4,4187473"^ 



July 

0-15 

8.585241“^ 

-1.439371“^ 

2.255601“^ 

-3.082077"^ 

2.079833"® 

-5.087808"^ 


15-47 

4.409444"^ 

-8.501794"^ 

6.408721“^ 

-2.137128"^ 

3.280708"^ 

-1.912325"® 


47-51 

-5.381483"^ 

3 . 443774 “^ 

-7.221806“^ 

4.994872"® 




51-65 

1 8.595289"^ 

[ 

-5.584211"^ 

1.371591“^ 

-1.5063“® 

6.229936"® 






for the CO 2 P20 line for the January and July 30 N model 
atmospheres. Their results for the optical depth 


, h » 

= J Pjj Q day as a function of terminus altitude for 


o -^2' 

vertical propagation from sea level are 


x(Jan) = 0.0475 [l-exp(-0, 705 h)] + 0.0877|^l-eKp(-1.15 h)j (1-58) 
and 

T(July) = 0.118[l-exp(-0.635 h)] + 0.429 [l-exp (-1.01 h)j (1-59) 

where h is the alttiude in kra. 

For LDV applications along atmospheric slant paths, 
one is interested in the two-way loss coefficient 

L(h2-hi) = 

~ cos 6 

5^ {■^BjO '*’2> - HjO'**!')] ■ 

These expressions can be readily evaluated using the relations 
given previously. 

Aerosol extinction is primarily due to scattering at 
wavelengths shorter than 6y. At longer wavelengths (at 
least in the continental aerosol) aerosol absorption is 
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comparable to scattering losses- Figure 1-15 shows the 

aerosol attenuation coefficients (optical depth per km 

16 

of path length) presented by McClatchey and Selby . These 
attentuation coefficients are for a ’’clear'* atmosphere at sea 
level with a particle size distribution comparable to that 
of Figure 1-13. They scale with particle density for other 
conditions. Table 1-4 gives particle densities as a function 
of altitude for two aerosol models, the "clear" atmosphere 
(Visibility 23 km) and a "hazy" atmosphere (Visibility 5 km) . 

In Table 1-5 we compare molecular and aerosol losses 
at O.Sy and 10. 6y for level paths at sea level and at 10 km 
altitude at different ranges and visibilities. These calcu- 
lations are for the July 30®N model atmosphere and the 5 km 
and 23 km visibility aerosol models. 



PARTICLE DENSITY 
N (PARTICLES PER cm3) 


Altitude 

(km) 

23 -km Visibility 
Clear 

5-km Visibility 
Hazy 

0 

2.828E+03 

1. 378E+04 

1 

1.244E+03 

5.030E+03 

2 

5,37 lE+02 

1. 844E+03 

3 

2,256E+02 

6,731E+02 

4 

1, 192E+02 

2*453E+02 

5 

8. 987E+01 

8.987E+01 

6 

6.337Et-0i 

6.337E+01 

7 

5.890Ef01 

5.890E+01 

8 

6,069E+01 

6.069E+01 

9 

5.8I8E+01 

5.818E+01 

10 

5.675E4*01 

5. 675E+01 

11 

5.317E+01 

5.317E+01 

12 

5. 58553^-01 

5.585E+01 

13 

5, 156E+01 

5* 156E+01 

14 

5.048E+01 

5,048E+01 

15 

4.744Ef01 

4.744E+01 

16 

4. 511E4-01 

4.511E+01 

17 

4.458E+0I 

4.458E+01 

18 

4.313E+01 

4.313E+01 

19 

3. 634E+01 

3. 634E+01 

20 

2. 667E+0i 

2.667E+01 

21 

1. 933E+01 

1.933E+01 

22 

1.455E+01 

1.455E+01 

23 

1. 113E+01 

1.113E+01 

24 

8.826E+00 

8. 826E+00 

25 

7.429E+00 

7.429E+00 

30 

2,238E+00 

2.238E+00 

35 

5.890E-01 

5.890E-01 

40 

1. 550E-01 

1.550E-01 

45 

4.082E-02 

4.082E-02 

50 

1.078E-02 

J.078E-02 

70 

5.550E-05 

5.550E-05 

100 

1.969E-08 

1.969E-08 








Table 1-5: 

Atmospheric 

Altitudes, 

Attenuation at Different Visbilities 
Ranges and Wavelengths. 

Altitude 

Visibility 

Range, 

2-way Attentuation in db 


(0.5y) 

km 

X=0.5y X=10.6y 

A* A* 

Aerosol Aerosols CO„ H„0 


Sea Level 


6.76 

20.3 

67.8 

204 

1.48 
4.42 
14. / 
44.2 
147 


0.62 

0.75 

1.84 

2.23 

6.16 

7.46 

18.5 

22.3 

0.13 

0.75 

0.40 

2.23 

1.34 

7.46 

4.02 ! 

22.3 

13.4 

74.6 
1 47 



.01 

4,27 

.03 

14.2 


For particle density of Table 1-4 


July, 30 N model atmosphere 






3.4 Frequency Properties of the Signal 


The range of beat frequencies Af = corresponding 
to wind velocities from 0 to 20 m/sec will be (0-80 MHz) 
cos$ at O.Sy and (0-4 MHz) cos(J» at lOy, where $ is the 
angle between the line of sight and the velocity vector. 

The frequency spread due to turbulent eddy velocities within 
the heterodyne volume will be independent of with a half- 
power width corresponding to the rms eddy velocity - 

Frequency broadening due to atmospheric turbulence 
in the optical path, calculated to be on the order of kilo- 
hertz in the preceding section, can therefore be neglected 
for any reasonable velocity resolution requirement. Broaden 
ing due to fluctuations in the laser frequency is larger (of 
the order of 20 to 100 khz for a non-stablized single mode 
laser) and could introduce errors of a few cm/sec into the 
velocity measurement. 

In situations where the total molecular scattering 
considerably exceeds the particulate scattering (for 
example, in maritime air) , interference by the thermal 
molecular motion may occur. About 5% of the molecular 
scattered signal will be present in the frequency range 
corresponding to velocities from 0 to 20 m/sec. If this 
component is comparable to the Mie scattering, it will be 
necessary to provide some sort of discrimination technique 
to separate the molecular and Mie compenents. 







The particle size distribution for maritime aerosols 

-I peaks at about l]i radius. Palmer and Zdunkowski^^ have 

estiraatedthe particle density necessary for the particulate 

backscattering from Ip particles at a wavelength of 0,7u to 

equal the molecular Rayleigh component. Scaling their 

results to sea level, we find that a particle concentration 
—3 

of 10 cm would produce backscattering equal to molecular 
scattering. Measured values of maritime aerosols are in the 

_3 

1 cm range. Integrated over frequencies corresponding to 
0 to 20 m/sec, the Rayleigh and Mie contributions to the 
power density spectrum are thus expected to be comparable. 
There are two significant differences in the signals, however. 
First, the power density spectrum due to molecular scattering 
will be essentially flat in the frequency range of interest 
so that scanning the spectrum would provide discrimination. 

Second, the angular distribution of wind velocities 
peaks in the direction of the mean velocity, whereas the 
thermal molecular velocities are essentially independent of 
the viewing angle. Thus, spatial scanning would also pro- 
vide a means of discrimination for systems operating in 
the visible. Molecular scattering are very weak for in- 
frared systems. 
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3.5 Choice of Wavelength 

18 

Bullrich has shown that the effective mean value 

of na for atmospheric aerosols in continental air varies 

inversely as the wavelength for X between 0.1 and 20y 

(see also Figure 1-15) . This conclusion is essentially 

dependent on the existence of an inverse cube dependence 

on the particle radius for dn/dlogr (see Figure 1-14) . 

For this type of particle distribution, the signal-to- 

noise ratio per joule output for an optimally designed 

system should exhibit a linear dependence on wavelength 

2 4 

in the near field and a X ’ dependence in the far field 
(assuming that comparable quantum efficiencies can be 
obtained at any wavelength) . In terms of the net laser 
output energy E (joules ) , the maximum signal-to-noise 
ratios are given by the expression . 


(S/N) = 1.25x10^^ n E C X 

ni3.x 


for L < L 

n 


= 4.2x10^® n E ? (6L/L) 


,12/5/^12/5 


^11/5 

Lt 


for L > . 

Here 5 is the wavelength independent scattering function 

na X and is the near field range. 

Thus, the figure of merit to consider in the selection 

12/5 

of a CW laser is nEX at short ranges and nEX at long 


ranges (with, of course, the requirement tiiat the detector 
be able to respond to the megahertz beat frequencies) . Since 




lr“i 




3. 6 Requirements on Detector ^Frequency Response 

In Table 1-6 we have tabulated the maximum detector 
time constants that can be tolerated for the detection of 
radial velocities in the 0 to 20 meter/second range for 
visible and IR wavelengths and for both a stationary plat- 
form and a platform moving at a velocity of 300 meters/ 
second. 


Table 1-6; Time Constants Required to Detect 
Atmospheric Motion in the Range 
0 to 20 m/sec (t = l/2frf) 


A 

0.5y 

lOu 


stationary 
Observer 
2 ns 
40 ns 


Observer Moving 
at 300 m/sec 
0.15 ns 
3 ns 


These rather stringent requirements for very fast 
detectors can be softened considerably (to time constants 
of the order l/2Tr Af = X/4 ttAv, where Av is the velocity 
resolution desired) if the frequency of the local oscillator 
can be suitably controlled. For Av of the order of 2 m/sec / 
this would permit use of detectors having time constants 
in the more easily obtainable range of 20 nanoseconds in the 
visible and 400 nanoseconds in the infrared. Separate tuning 
of the local oscillator can be obtained by using separate, 
frequency stabilized, and tunable lasers for the transmitter 
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and the local oscillator. Alternately, single or multiple 
reflection of the reference beam from a moving mirror will 
permit sweeping of the local oscillator frequency. This 
latter technique is most useful when sinusoidal repetitive 
sweeping of the local oscillator frequency is acceptable. 
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In a typical coaxial heterodyne system, both the 
transmitted and received beams suffer a 50% loss at the beam 
splitter. In a side-by-side system, these losses can be 
reduced by using beam splitters with transmission > 50%. 

The following calculations will thus be somewhat conservative 
for side-by-side systems. Some sort of modulating chopper 
is assumed to achieve the narrow band detection corresponding 
to integration times on the order of seconds. An efficiency 
of 25% in the remaining optical and electronic components 
would result in an overall optical-electronic efficiency 
in the neighborhood of 3%. 

Assuming 3% optical efficiency, 20% quantum efficiency 

and a range resolution requirement of 10%, the signal-to- 

noise per joule of transmitted energy is presented in 

Figure I- 16 as a function of range in maritime and continental 

atmospheres at wavelengths of 0.5 and lOp and for turbulence 

“8 “X/3 

levels of 10 and 10 cm ' . Atmospheric attenuation has 
been neglected for the sake of clarity. 

In order to be able to measure the mean wind velocity 
to within 10%, it is estimated that the power S/N should 
be at least 300. Using a 10 watt CW Argon laser operating 
at 0.48y, the daytime detection range is thus seen to be 
on the order of 100 meters in a maritime atmosphere and a fev/ 
hundred meters in continental air, assuming an integration 


S/N PER JOULE 








time of one second. A 10 watt CO 2 laser (10. 6y) would 
extend the effective daytime measurement range to something 
like 10 kra. Coupled with the fact that 10 watts is currently 
the upper limit of available CW lasers in the visible, 
with efficiency (ratio of power out to power in) consider- 
ably less than 0.01%, whereas CO 2 lasers provide up to one 
kilowatt CW power with efficiency better than 10%, it would 
appear to be strongly desirable to go to 10. 6 y. 

Disadvantages of operating at 10. 6 p include the 
necessity of cooling the detector (best detector available 
with adequate frequency response is Ge;Hg, which must be 
cooled to less than 40® K) , a requirement for CO 2 # ^ 2 , and 
He flow for the laser (although sealed units are now be- 
ginning to appear on the market) and difficulty in aligning 
the optical system outside of the visible region. 

In most applications, the advantages of a CO 2 system 
will outweigh the disadvantages. For long-term untended 
operation, however, the requirement for liquid He (or H 2 ) 
may preclude the long wavelength system. 

To evaluate the effect of a physical limitation on 
the size of the optical system, which may be the case in 
some practical applications, the two upper curves in Figure 
1-16 illustrate the degradation in S/N as a function of range 
if the aperture diameter is limited to 30 cm. For daytime 
turbulence, the effect is negligible; for nighttime turbu- 
lence levels, the S/N may be reduced by as much as an order 
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of magnitude from the S/N achievable with no size limita- 
tions on the optics. The effect is more pronounced at 
altitude where lower turbulence levels permit the use of 
large apertures. Figure 1-17 illustrates the S/N per joule 
as a function of range at an altitude of 10 km. If one 
assumes an upper limit of 30 cm for the aperture diameter 
(which would be a reasonable upper limit for aircraft) , the 
S/N advantage of the longer wavelength is negated by the 
size restriction. 


In Figure 1-18 two of the curves from Figure 1-16 
are reproduced and the effects of atmopsheric attenuation 
are included. From the figure it can be seen that atmo- 
spheric attenuation becomes important at ranges where the 
S/N is becoming marginal, putting rather sharp limits on 
the detection range possible. 


0 ® 

RANGE (METERS 


Figure I~l7 Signal-to-Noise per Joule vs. 

for Horizontal Path at 10 km 


VALUE FOR 30 cm. 
APERTURE DIAMETER 






Figure 1-18- Signal-to-Koise per Joule vs. Range 
(Atmospheric attenuation included) - 



(SIGNAL /NOISE PER JOULE) 



'a.;' 



3. 8 CW versus Pulsed Doppler 'Lidars 

The theoretical limits for the S/N values that can 
be achieved in an optimally designed heterodyne system are 
the same for both pulsed and CW operation so long as dif- 
fraction limited antenna patterns can be achieved in both 
modes. Thus, the detectivity limit for ideal systems 
is dependent only on the average output power level. Large 
differences can be expected, however, in practical imple- 
mentation. The major advantages of a pulsed system are 
that a single coaxial system can be used at all ranges and 
that range resolution is obtained automatically with time, 
whereas in a CW system at long ranges, separate transmitter 
and receivers are required and a precision mechanical adjust- 
ment of the scattering angle is needed to change ranges. 

Also, a pulsed system is much less sensitive to the detailed 
structure of the antenna pattern, whereas a CW system can 
suffer greatly from the existence of side lobes. Nevertheless, 
if greater mean, single mode, diffraction limited output 
powers can be obtained in CW operation, it is clear that this 
mode would have the greatest inherent sensitivity. Thus, 
a CW system coupled with amplitude or frequency modulation 
or coding techniques to obtain the ranging advantages of a 
short pulse system may be the more appropriate for practical 
applications . 
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3-9 A NOTE ON PREDOMINANT NOISE IN THE HETERODYNE SIGNAL 

In the preceding discussion, the limiting noise in the 

heterodyne signal was assumed to be the photon noise in the 

reference beam. The conditions under which this assumption 

is justified are as follows - 

In most heterodyne systems, the reference beam is 

attenuated to avoid detector saturation and to minimize the 

effects of any intensity fluctuations in the reference beam 

which might occur at the heterodyne frequency. If, after 

this attenuation, the signal is still limited by shot noise 

in the reference beam, the S/N will be independent of the 

actual value used for the attenuation. 

The noise terms which must be evaluated include the 

• — X 

photon noise in the scattered signal, given by n N T no 

fluctuations in the laser output at signal frequencies, given 

by a n N T 6 , where a is the attenuation of the reference 

beam and 6 the fractional fluctuation per unit bandwidth, 

detector noise given by 2 i^ t and background noise. The 

worst background condition will be that for sunlit clouds, 
n Nj^ AX T 

given by ^ A^, where N,^ is the spectral radiance 

of the clouds , the acceptance angle and A the •'■ollecting 

area of the receiver optics. For a diffraction limited 
system the product of the solid angle of acceptance and the 
collecting area is equal to X . 
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A comparison of these terms and the implied restrictions 
on the attenuation coefficient a are given in Table 1-7. For 
strong signals, i.e., when the scattered intensity exceeds 
the backgrotmd level and the detector NEP, a is determined 
as follows: 


0 


6 n N T 


>> a >> na 


For the example in Table I- , 


^ » a >> 6.5x10"^^ 


This is equivalent to stating that the S/N given in 
Figure 15 cannot be achieved if fluctuations in the laser 
output approach 1% in the frequency range of interest. 

A comparison of background and detector noise terms 
in Table 1-7 indicates that detector noise will exceed back- 
ground noise (for the example chosen ) , even for the O.lu 
spectral bandwidth assumed in the example. Background levels 
can be further reduced by orders of magnitude by means of 
intereference filters or monochromators if necessary. Thus, 
for weak signals, a is determined by detector noise and 
laser output fluctuations: 


6^ n N T nN 


147 


Table 1-7: Conditions .for Photon Noise in 

Reference Beam to Predominate 


Source 


Noise Term Necessary Conditions 


General Case 


Reference Beam a n N x 


Example 


Scattered Beam n N x na ^ 


a >> na 1 - 


a >> 6.5>«10 


-13 


• 2 1 6^10 
Laser Output (a 6 n N x) a << —= ; — a << 

Ml I * 


Fluctuation 


6 n N X 


-17 


2 


Cloud Radiance X 


a >> 


N hv 


9 —14 

a » 6.5X10 


Detector 


2id " 


a >> 


2i 


nN 


a >> 8X10 


-13 


For following conditions: 


X = 0.5U 


-1 2 

N-i = 10 w/cm ster y 


n = 0.2 


AX = O.ly 


N 


= = ^.SxlO^"^ ph/sec i^ = 


4 -14 

6x10 elec/sec (10 amps) 


X = 1 sec 


na 


1.3x10 ^ m"^ 
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or 


~ >> a n N T >> 2i- T 


6 


For = 6x10^ 
should be not more 


-14 

elec/sec (10 amps) 
than a few tenths of a 


and T = 
percent 


1 sec. 
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APPENDIX II 

Doppler Offset 

One of the principal problems with the current 
airborne system is a ringing effect that allows the trans- 
mitter to continue to transmit small amounts of energy 
after its nominal cutoff time. 

The Doppler frequency offset in a ground-based 
system is small and some method for isolating the receiver 

from scattered radiation from the transmitter is required. 

6 10 

For a range of 10 km the round-trip time is 2x10 /3x10 
= 70y sec. A mechanical shutter can easily operate in this 
time, i.e., a disk rotating at 600 rps (36,000 rpm ) and 
having diameter of 4 cm has peripheral speed 271(600) *2 = 
7539 cm/sec. In 50y sec a slot at the edge moves 
7539x50x10 cm = 0.39 cm = 3.7 mm. Thus, by focusing the 
beam to a small spot on such a switch and synchronizing the 
transmitter, the transmitter feed-through can be strongly 
blocked. Two such switches could be used to eliminate 
virtually all feed through. At a range of 2 km the round- 
trip time is 14iJ sec (the time after start of pulse that 
signal begins to return from 2 km) . Here we would need a 
switch that operates in less than 14y sec, i.e., a slot 
width of 0.7 mm. 
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APPENDIX III 


In this section we shall derive the dispersion 
relation for the Helmholtz waves in a confined system of 
two steadily flowing layers of different densities. Al- 
though it is possible to derive the formula for a general 
three-layer model with the middle layer of higher stability, 
it is more straight-forward to deal with the asymptotic 
case when the depth of the middle layer approaches zero 
and both the top and bottom layers are of constant densities. 
Because of the constant densities, if the flow is initially 
irrotation 'i, it will remain so. Velocity potentials 

and exist in both the top and bottom layer. Figure A 
shows the coordinates and the wave form of the interface. 


The equations of motion in the small perturbation 
limit are: 


P« u 


Po « 


Pi U 


Pi u 


9^0 _ 

3Po 

3x 

3x 

3w _ 

9Po 

3x 

3z 


3P 


9x 

3w, 


3x 


Po 9 


3Pl 

^ - Pi 5 


(III-l) 


(III-2) 


(III-3) 


(III-4) 


'1 '' 3x 

where subscript o, 1 denote the corresponding variables in the top 







and bottom layer and u, w, p, p' are the perturbed variables. 

Since u^ = Vti)^ and the pressure can be expressed in 

terms of cji^ and (j>^ by integrating (III-l) and (III-4) with respect 
to X and z correspondingly. Hence, we obtain: 

3<l>o 

Po = - Po ^ 

Pi “ ^ Pi ^ "3F 

Assuming a standing wave 
at the interface so that 

n = a sin(kx) 


(III-5) 

pj[ gz . (III-6) 

of the unknown wave number k exists 
the displacement is given by 

(III-7) 


we apply the boundary conditions: 


kinematic condition: 




U 


9IL 

3x 


n 


(III-8) 


dynamic condition: the pressure in each layer given 

by (III-5) and (III-6) should agree at the interface. 


Further boundary conditions on the top and bottom walls are 
that normal velocity should vanish, i.e.. 


3z 


h 


0 


3 z 


-h 


0 


(III-9) 


(III-IO) 
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The general solutions to V <j> = 0 are 

(j)^ = coskx ^ A exp(kz) + B exp(-kz) j 

<(>j^ = coskx [c exp(kz) + D exp(-kz) ) 

Applying the boundary conditions (III-5) , (III-6) and 
we can obtain the following relations: 

Uk(A+B) = pj^Uk(C+D) - pj^ga 

A-B = C-D = Ua 

Equations (1II-9) and (III-IO) require 
A exp(kh) “ B exp(-kh) = 0 
c exp(-kh) - D exp(kh) = 0 

Eliminating A, B, C, D from Equations (III-13) through 
we obtain the relation determining the k. 

k - —= tanh(kh) = 0 

u"(PltPo) 


(III-ll) 

(III-12) 
(III-8) , 

(III-13) 

(III-14) 

(III-15) 
(III-16) 
(III-16) , 

{III-17) 


g(pi"po) _i 

Substituting — .1027 km and h = 10 km, the 

° <Pl+Po> 

root for Equation (17) can be found to be .025 km . This 
agrees well with the numerical result .02576 km ^ obtained 
from the present code solving a three^layer model. This 
model is composed of a thin middle layer with a depth one 
thousandth of that of its neighboring layers. The 
middle layer stability is 4.52 km"^ which is 1369 times 
higher than that of both the top and bottom layer. 


I 

j 


1 
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